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We suggest a new method of calculations for a clean Fermi gas with a repulsion in any dimension. 
This method is based on writing equations for quasiclassical Green functions and reducing them to 
equations for collective spin and charge excitations. The spin excitations interact with each other 
and this leads to non-trivial physics. Writing the solution of the equations and the partition function 
in terms of a functional integral over supervectors and averaging over fluctuating fields we come to 
an effective field theory describing the spin excitations. In some respects, the theory is similar to 
bosonization but also includes the "ghost" excitations which prevents overcounting of the degrees of 
freedom. Expansion in the interaction reveals logarithmic in temperature corrections. This enables 
us to suggest a renormalization group scheme and derive renormalization group equations. Solving 
these equations and using their solutions for calculating thermodynamic quantities we obtain explicit 
expression for the specific heat containing only an efi^ective amplitude of the backward scattering. 
This amplitude has a complicated dependence on the logarithm of temperature, which leads to a 
non-trivial temperature dependence of the specific heat. 

PACS numbers: 71.10.Ay,71.10.Pm 



I. INTRODUCTION 

Landau theory of the Fermi liquid (FL) suggested 50 
years ago^ is the basis for description of normal metals. 
Roughly speaking, the main statement of the FL theory 
is that the low energy behavior of interacting fermions 
is similar to that for the ideal Fermi gas. This theory 
explains very successfully properties of a large number 
of metals and He'^. It is quite common to discuss ex- 
perimental systems just forgetting about the interaction 
and using phenomenological parameters like, e.g. effec- 
tive mass, density of states at the Fermi surface, etc., 
instead. 

Yet, a recent progress in study of unconventional 
metals like high temperature superconductors, heavy- 
fermion materials has revealed considerable deviations of 
their properties from those predicted by the FL theory 
(for a recent review see, e.g. Refs. 2,3. As a result, quite 
a few theoretical works have appeared recently where the 
validity of the Landau FL theory was discussed^"^. 

At first glance, being phenomenological from the be- 
ginning, the Landau theory has been later confirmed by 
analyzing diagrammatic expansions^" (see, also Refs. 11- 
14) and looks very well established. However, a very 
strong assumption was used in this discussion, namely, 
that one could single out a singular particle-hole chan- 
nel and sum proper ladder diagrams. Two particle ir- 
reducible vertices entering the ladder diagrams should 
remain finite and analytic in the limit of small momenta 
and frequencies. Of course, this is not always true and 
under certain conditions the system may become super- 
conductor, antiferromagnct, etc. In many cases, the fail- 
ure of the Fermi liquid description can be checked by 
a more careful consideration of the perturbation theory 



and, e.g., the existence of the superconducting transition 
can be established in this way^^"^'^. 

Nevertheless, it is believed that the system of fermions 
with a repulsion should behave like a Fermi liquid pro- 
vided the dimensionality d > 1 and there are no van-Hove 
type singularities on the Fermi surface. Naturally, the 
similarity between the Fermi liquid and ideal Fermi gas 
cannot be exact and, clearly, there are corrections at fi- 
nite temperatures, finite frequencies or momenta. These 
corrections become especially interesting when they are 
non-analytic functions of the values of the temperatures, 
frequencies or momenta. 

For the ideal Fermi gas, such quantities as C (T) /T 
and x{T), where C (T) is the specific heat and x{T) 
is the spin susceptibility, can be represented in a form 
of asymptotic series in T'^/e'^ {ep is the Fermi energy). 
Fermion-fermion interactions lead to additional contribu- 
tions to these quantities that are not necessarily analytic 
in . It is well established that in D = 3 next-to- leading 
term in C (T) /T is T^lnT, see Refs. 15-19. It was 
claimed in Ref. 20 that the non-uniform spin susceptibil- 
ityi x{Q)^ depends on the momentum Q as Q^l'n.Q. In 
2D, non-analytical corrections to C (T) /T and x {Q: T) 
found so far scale as T, see Refs. 21-24, and max {Q, T}, 
see Refs. 20,22,25-28, respectively. 

The existence of the non-analytical corrections to the 
physical quantities is not accidental. In fact, all of the 
singular corrections to the thermodynamic quantities can 
be understood in terms of the dynamics of the low lying 
collective excitations, see e.g. Ref. 24,28. A detailed 
analysis of non-analytic corrections to the specific heat 
of a three-dimensional Fermi liquid is given recently in 
Ref. 18. It was shown in Ref. 28 that all the other contrib- 
tutions contain integrations over the entire Fermi surface 



2 



arc regular in T^, unlike the contributions of the collec- 
tive modes which contain 2kp scattering. 

Explicit calculations for systems such low lying modes 

arc not simple even in the lowest orders of the perturba- 
tion theory. This situation is analogous to that in theory 
of disordered metals, where the low energy behavior of 
the system is governed by the multiple interference of the 
electron waves scattered by impurities. In the diagram- 
matic language, this effect can be expressed in terms of 
an interaction between electrons and diffusion modes (so 
called cooperons and diffusons^^'^'^). Calculations in high 
orders in the diffusion modes (weak localization correc- 
tions) using the diagrammatic expansions are also quite 
involved. 

However, another approach has been developed in the 
theory of disordered metals based on integrating out elec- 
tron degrees of freedom and deriving an effective La- 
grangian describing the diffusion modes. This reduction 
simplifies calculations because only low lying excitations 
are left in the theory. The Lagrangian has the form of a 
so called cr-model, first introduced in the theory of disor- 
dered metals in Ref. 31 using the replica trick. Another, 
supermatrix form of the a-model, is based on a supervec- 
tor representation'^^ of Green functions, and this method 
has found numerous applications (for a review, see, e.g., 
Ref. 33). 

One can see a certain analogy between calculating the 
non-analytical corrections for the Fermi gas with inter- 
action and the weak localization corrections in theory of 
disordered metals. Following this analogy it seems quite 
natural to try to develop a scheme that would allow us 
to reduce the initial model of the interacting Fermi gas 
to a model describing only low lying excitations. Then, 
we would have in the theory not the initial fermions but 
the collective excitations like, the zero sound, may be, 
weakly interacting with each other. Apparently, the lat- 
ter are bosons and we would have, as a result, a system 
of bosons instead of the initial fermionic system. Reduc- 
ing fermion models to boson ones is usually referred to 
as hosonization and we will loosely use this word in the 
subsequent discussion, though the final theory we derive 
will be necessarily supersymmetric, and contain also the 
fermionic degrees of freedom - "ghosts" . 

The notion of the supersymmetry appears quite nat- 
urally already in the consideration of the leading non- 
analytic corrections to, say, specific heat, C . Namely, for 
the spinless electrons, the ring-diagrams correction to C 
is of the form 



ization operator acts on any function 6(n, q) as 



SC = -T 
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where the unit vector n characterizes the position of the 

momentum on the Fermi surface, and a;„ = 2TmT is a 
bosonic Matsubara frequency. The quasiparticle polar- 
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with vp being the Fermi velocity. The operator F is de- 
fined as 



F6 



(ni, k) = y du2¥ (nin2) 6(n2, k) 



and F (111112) is the Fermi-liquid function describing the 
interaction between the quasiparticles moving in direc- 
tions 111 and 112. We can imply the proper normalization 
of the solid angle, i.e. 
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drv = 1, 



(1.2) 



when integrating over the momentum directions 
Equation (1.1) can be re- written identically as 
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[In all the consideration we will not write factors of 
volume, whenever they are self-evident for an educated 
reader] . The form of this expression is quite instructive. 
The first term, D.p is nothing but the contribution of the 
non- interacting bosons, whose spectrum is determined by 
the kinetic equation in Landau theory. However, those 
bosons are made out of electrons which are already in- 
cluded in the leading term of the specific heat. That 
is why the second term, fig, simply subtracts the con- 
tribution of the electron-hole pairs in the absence of the 
interactions to avoid a double counting. Since the contri- 
bution of the second term is opposite to the contribution 
of the physical bosons, it will be natural to treat them as 
pseudofermions, or "ghosts" and include them into the 
field theory description on equal footing with the bosonic 
field'. 

However, one cannot proceed in a direct analogy with 
the supersymmetry method of Ref. 33 because the lat- 
ter is essentially based on the use of a sufficiently strong 
disorder leading to a diffusion motion at not very long 
distances. Fortunately, the method can be generalized 



' The partition of the low-energy excitations in the Fermi liquid 

in terms of physical bosons and artificial ghosts (to avoid over- 
counting of the degrees of freedom) was suggested in Ref. 24, 
however, the consideration there was limited to the theory of 
non-interacting bosons only. 
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to the clean case by writing equations for quasiclassi- 
cal Green functions and representing their solution in 
terms of functional integrals. Equations for the quasi- 
classical Green functions for the disorder problems were 
introduced in Ref. 34. The authors of Ref. 34 suggested 
writing their solutions from the condition for a minimum 
of functional having a form of a ballistic cr-modcl. This 
could be written in a form of a functional integral pro- 
vided this integral could be calculated by the saddle point 
method. At the same time, it was not clear why the sad- 
dle point approximation could work well for the clean of 
weakly disordered case. 

Later it was realized^^ that the solution of the qua- 
siclassical equations can be exactly written in terms of 
a functional integral with a Lagrangian having the form 
of the ballistic a-model. Within such an approach one 
reduces the initial electron model to a model describing 
"collective excitations" and therefore it is relevant to call 
it loosely bosonization. A smooth potential could be con- 
sidered in this approach but no interaction was included. 

In the present paper we develop a new method for 
studying clean fermion systems with a fermion-fermion 
repulsive interaction. It is based on decoupling the in- 
teraction by a proper Hubbard-Stratonovich transforma- 
tion. Both the forward and backward scattering are 
taken into account and therefore the slow decoupling field 
$ (r,T) has a spin structure (where r is the coordinate 
and T is the imaginary time). After the decoupling we 
derive equations for the quasiclassical Green functions. 
In order to solve the equations we use a new trick based 
on a slow dependence of the field $ (r,T). As a result, we 
obtain linear non-homogeneous equations for spin and 
charge excitations. The spin excitations are most impor- 
tant and we represent the solution of the corresponding 
equation in terms of an integral over supervectors. Aver- 
aging over the field we obtain an effective theory with 
a V"'' + V-"^ interaction. This theory describes collective 
bose excitations. It is important to emphasize that the 
ip^ and tp^ terms arise due to spin-spin interactions. For 
fermion models containing only a density-density inter- 
action, those terms vanish and the theory becomes free. 

Making expansions in the tp^ + ip^ interaction we found 
that the theory is logarithmic in any dimension, which 
allowed us to use a renormalization group approach to 
sum up the parquet series. Writing and solving the RG 
equations we are able to express the non-analytical con- 
tribution to the specific heat with a logarithmic accuracy. 

It is relevant to mention that the method we develop 
now is completely different from the u-model approach 
for disordered systems with interaction^^ (see also a re- 
cent siipersymmetric formulation''^). In the present ap- 
proach the collective excitations are described by a La- 
grangian with the ip^ + ip'^ interaction and not by a a- 
model. 

Attempts to bosonize fermionic models in the dimen- 
sionality d > 1 have been undertaken in the past start- 
ing from the work''^. In this first work the one dimen- 
sional bosonization was directly extended to higher di- 



mensions. This idea was hirther developed more recently 
in a number of publications^''^''^'^^^®. Our approach is 
completely different and more general. The high dimen- 
sional bosonization developed previously can be applica- 
ble only for a long range interaction when the backward 
scattering is absent. In this case the fermion-fermion 
interaction is replaced by an interaction of the local den- 
sities (local in space and in the position on the Fermi sur- 
face). This means that effects related to electron spins 
are beyond the possibility of that method. In contrast, 
the backward scattering is included in our approach and 
plays a very important role. 

The article is organized as follows: 

In Section II, we formulate the model and single 
out slow pairs in the interaction term. We perform 
a Hubbard-Stratonovich transformation and reduce the 
model with the interaction to a model with slowly fluctu- 
ating fields. Then we derive quasiclassical equations for 
electron Green functions. 

In Section III, we introduce an eikonal type method 
(a.k.a. Schwinger Ansatz) for solving the equations. As 
a result, we obtain equations for effective charges and 
spins in the presence of fluctuating fields. We express 
the partition function in terms of the solutions of these 
equations and calculate it neglecting the interaction be- 
tween the excitations. 

In Section IV, we write the solutions of the equations 
for the collective modes and the partition function in 
terms of functional integrals over supervectors. We av- 
erage over the fluctuating fields and derive an effective 
field theory containing interaction terms. 

Section V contains an explanation how one obtains log- 
arithmic contributions. Then, we develop a RG scheme 
integrating over fast variables and writing renormalized 
coupling constants. 

In Section VI, we derive renormalization group equa- 
tions and find their solutions. 

Section VII is devoted to calculation of the thermody- 
namic potential and specific heat using the solutions of 
the RG equations. Explicit formulae are obtained in two, 
three and, separately, in one dimensions. 

Our findings are discussed in Section VIII. 



II. THE MODEL AND BASIC EQUATIONS. 

A. Formulation of the model and singling out slow 
modes 

We start with formulating the model we would like to 
investigate. This is the most general model for fermions 
with a short range interaction " in an arbitrary dimen- 



" Inclusion of the long-range Coulomb interaction is straightfor- 
ward and does not lead to any consequences relevant for our 
study 
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sion d. The Fermi surface is assumed to have no singu- 
larities. In order to avoid unnecessary trivial generality 
we consider the Fermi surface to be a just d -dimensional 
sphere. 

It will be convenient for us to express the physical 
quantities in terms of a functional integral over anticom- 
muting variables x with an Euclidian action S. In this 
formulation, one can use imaginary time r in the interval 
< r < 1/T, where T is the temperature and write the 
temperature Green function G^a' {x,x'), as follows 

G,,,, {x, x') = J xa {x) xl- {x') exp {-S) DxDx* , 

(2.1) 

where 

x={r,T), jdx-=jd'^vjj dT- (2.2) 

r is the spatial coordinate, and a labels the spin. 

The partition function Z entering Eq. (2.1) has the 
form 

Z = j e^p{-S)DxDx*. (2.3) 
The action S in Eq. (2.1) can be written as 

S = j Codx + Sint, (2.4) 
where the term Cq, 

■Co = E (^) (- ^ - ^o) Xa (x), (2.5) 



Ho = sf 

2m 



(2.6) 



stands for the Lagrangian density of free fermions {ep is 
the Fermi energy, m is the mass and p is the momentum 
operator) and Sint describes the fermion-fermion inter- 
action, 



(T,<t' 

X [xl (x) xl' {x') Xa' {x') Xa {x)] 



(2.7) 



where v{x — x') = U [r — r') 5 {t — t') and [/ (r — r') is 
the potential of the interaction. 

The field variable x must be antiperiodic in r with the 
period 1/T 

x(r,T) = -x(r,r + l/T). (2.8) 

The thermodynamic potential il can be written as 

n = -T\nZ. (2.9) 

The functional integrals over x with the Lagrangian 
C, Eqs. (2.4)-(2.7), are too complicated to be calculated 



exactly and making controllable approximations is in- 
evitable. For performing further formal manipulations 
we restrict ourselves with the case of a weak interaction, 
and discuss the changes of the theory for stronger inter- 
actions in the end of Sec. IIIB. As we have mentioned 
in the Introduction, the most interesting contributions 
come from the interaction of the fermions with low lying 
collective excitations and we would like to concentrate 
on such contributions. Thus, wc will try to simplify the 
interaction term Sint to display these collective modes 
explicitly. 

This can be achieved by singling out in the interaction 
term Sint pairs of the variables x slowly varying in space. 
Using the Fourier representation we write the effective 
interaction Sint containing the slow pairs as 



Sint^Sint = \Y. j dPidP2dK (2.10) 

cr,cr' 

X {y (k) xl (Pi) Xa {Pi + K) xl' {P2) Xa' {P2 - K) 

- V (P12) xl {Pi) Xa' {Pi + K) xl, {P2 + K) Xa {P2) }, 
where 



F(p) = j dre-'P'"t/(r), 

andpi2 = Pi-P2- InEq. (2.10), Pi = (pi,e„J , where pi 
is the momentum and = ttT {2ni + 1) arc Matsubara 
fermionic frequencies {i = 1,2). Short hand notation K 
reads K = (k, a;„), where a;„ = 27rTn are Matsubara 
bosonic frequencies. 

The symbol of the integration J dPi in Eq. (2.10) reads 

as 



jdPi{...)=TY,j 



d'^p 



(27r) 



(2.11) 



whereas the symbol / dK has the meaning 

/.K(...,=Ti://(k)£^^)^ (2.12) 



(2»)" 

In Eq. (2.12) we define the cut-off function 
/ (k) = /o (fcro) , fc = |k|. 



(2.13) 



The function /o (t) has the following asymptotics: 
/o (t) = 1 at f = and / (t) ^ at t ^ 00. This 
function is written in order to cut large momenta k. The 
parameter ro is the minimal length in the theory and we 
assume that ro much larger than the Fermi wavelength 
Xp = i/pF- In other words, the momenta k are cut by 
the maximal momentum kc = r^^ <^ pp, and the parti- 
tion (2.13) is not threatened by double counting, see also 
the discussion in the end of this subsection. 
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Equation (2.10) permits the further simphfications for 
the short range interaction potential. In the first term 
one can neglect the dependence on the transmitted mo- 
mentum k, 



F(k«pjr) = y2. 



(2.14a) 



In the second term one notices that momenta pi2 are 
close to the Fermi surface, so one can write"' 

V (pi2) = Vi (pTp^) ; Vi {9) = V (2pF sin ■ 

(2.14b) 

In Eq. (2.10) we recombine the terms with the help of 
the identity 

where cr = iox-'^yi'^z)-, and cXi, i = x^y,z are the 
Pauli matrices. Utilizing definitions (2.14), we re-write 
Eq. (2.10) as 



Sint = 2 y dPidP2dK\p {Pi,K) Vs {612) P {P2,K) 
+ Si{PuK)Vt{9i2)Si{P2,K) 



i=x,y,z 



(2.15) 



where, as before, 9i2 = PiP2- and the definitions 

p{p,K) = j:x: {p+§)x. (p-f) 

S, (P, i^) = ^ ar^xl, (P + f ) X.. (P - f 

(71 ,(72 ^ ' ^ 

(2.16) 

are introduced. The functions Vg (9) and Vi (9) are known 

as amplitudes of the singlet and triplet scattering, respec- 
tively. They are related to the amplitudes Vi and V2 as 

Vs{9) = V2-lvm, Vt{9) = -lv^{e). (2.17) 

In what follows we assume the operator Vs{9) 



Vs,tb (ni) = / dnaK.t {i^2) ^(na). (2.1^ 



to be positive definite and the operator Vt (9) to be nega- 
tive definite. According to Eq. (2.17), those assumptions 
imply the repulsive interaction. 



Equation (2.15) can be recast in a more transparent 
form. Performing Fourier transform over K in Eq. (2.15), 
we obtain 



i / dPidP2 1 dr f ^ dr 



X \p{Pv,v,t)Vs{9i2)p{P2;v,t) (2.19) 
+ Si{Pv,T,T)Vt{9i2)Si{P2;r,T) 



Wc borrow the onc-dirnensional notation of Ref. 47 for the cor- 
responding amplitudes. 



The entries in Eq. (2.19), 
p{P;r,T)=TYj ^e*--"V'/^(k)p(P,if); 

Sj{P;r,T)=TY J ^e*--"V'/^(k)5,(P,i^), 

(2.20) 

have the meaning of the smooth charge and spin den- 
sity accumulated in the phase space, and Eq. (2.19) is 
equivalent to the Landau description of the interacting 
quasiparticles. Appearance of the cutoff function /(k), 
Eqs. (2.12), (2.13), means that those densities may vary 
only with the spatial scale much larger than the Fermi 
wavelength Xp- 

Equations (2.15) - (2.20) constitutes the reduction of 
the original interaction to the interaction involving the 
soft electron-hole pair only. These are the only terms 
that may produce the non-analytic contributions to the 
observable quantities. In what follows we will manipulate 
with interaction (2.19) to obtain the low-energy theory 
in terms of the charge and the spin densities in the phase 
space. 

Closing this subsection, we discuss a very crucial is- 
sue that might start worrying an attentive reader at this 
point - what is the fate of the Cooper channel? In- 
deed, examination of the scattering processes induced by 
the Hamiltonian (2.10), see Fig. 1 shows that the vertex 
V3, describing the particle-particle interaction between 
the pairs x(Pi)x(— Pi) with x* iP2)x* {~P2) is missing. 
These are just terms that generate Cooperons for the 
system with the time reversal symmctry'^'^'''®. 

At first glance, we might miss important contributions 
because the Cooper channel generates logarithms in any 
dimension^^ and it looks as if we neglected them. How- 
ever, this is not so, although the reason for necessity to 
neglect the third vertex V3, Fig. Ic is rather non-trivial. 

As we will see, the most interesting contributions to 
interaction vertices and, finally to the thermodynamic 
quantities originate from the scattering on angle either 
or TT and not from an integral over the entire Fermi 
surface, so only such scattering amplitudes will be im- 
portant. 

Therefore, we must investigate the effect of the inter- 
actions in the Cooper channel on these particular ampli- 
tude, sec Figs. 2 Direct comparison of those con- 
tributions with the diagram generated by vertex Vi, see 
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Pl,fTl 



Pi + K, cri 
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Pi 



P\ 



Pi, (71 f Pi — K,(ji —Pi,ai —P2,cn 




Pl+K,a2\ P2,(J2 Pi + if, (72 ( P2+K,a2 




-P,(72 



-Pi 



-P',0-2 



FIG. 1: Vertices describing separation into soft modes, see 
Eqs. (2.10) or (2.19). Momenta k are small, k < kc <. pp. 
Diagrams a,b) correspond to the first and second terms in 
Eq. (2.19). The vortex c) was rightfully omitted in the low- 
energy effective theory, see text. 



Figs. 2b, 3b, shows that the perturbation theory in V\ 
generate in particular all the terms ''^ the Cooper chan- 
nel is responsible for. (These analogies in the two lowest 
orders of the logarithmic expansion are easily followed to 
all the higher orders). The only difference is the region 
of the integration over the intermediate momenta. In the 
diagrams of Figs. 2a, 3a, integrals over the intermediate 
momenta pi,P2 span all over the Fermi surface. On the 
other hand, (due to the condition |k| ~ 1/ro <C pf), the 
change in the momenta in Figs. 2b, 3b is very close to 
the backscattering. To be concrete, Fig. 3b) describes 
the region of momenta 



PPil 



P2P1I 



P2P'| < 



in Fig. 3a), where the maximal angle appearing in our 

theory is ^* ~ (pfT"o)~ • 

Taking into account the third vertex V3 in such region 
of the momentum space would mean double counting the 
contribution of this region that is most important for our 
consideration ^ . On the other hand, the region with the 



In fact, those terms will be properly accounted for in the renor- 
malization group treatment of Sees. V and VI. 
This double counting would not be important if all other regions 
of the Fermi surface contributed as well. In this case an error 
due to the double counting would be of order (roPF)^ '* <S 1- 
In dirty samples, the scattering on impurities leads to such an 
isotropization and therefore one should include in the effective 
Lagrangian the third term with the vertex V3 . The possibility to 
neglect it is specific for the clean systems considered here. The 
same approximation has been used in Rcf. 22 for computation of 
diagrams of the conventional perturbation theory. 




-P-K,a2 



-P-Ki 



FIG. 2: a) Lowest logarithmic contribution due to the Cooper 
channel, ai ^02- b) The same contribution obtained as an 
interaction of the spin modes generated by vertices Fig. lb, 
and fe -C pf- This diagram coincides with the renormalization 
of the quadratic part of the spin wave Lagrangian, see Fig. 9b) 
re-written in terms of electronic lines. 



scattering angle exceeding does not appear very inter- 
esting. One could integrate over this region from the be- 
ginning and this would simply renormalize the coupling 
constant Vt. As we work in the limit of a weak coupling, 
this renormalization cannot lead to any non-trivial effects 
and we assume in the subsequent discussion that it has 
already been performed. We will return to the discussion 
of this renormalization in Sec. VII C. 

One more uncertainty in the channel separation ap- 
pears when all four momenta in Sint of Eq. (2.7) are close 
to each other. In this limit the separation into slowly 
varying pairs in Eq. (2.10) is ambiguous. We resolve this 
unccertainty by attributing this region of the momenta 
to the vertex V\. As a consequence, the vertex V2 of 
Eq. (2.14a) should vanish for the momenta pi and P2, 
such that |pi — P2I — ^ Vf- However, the vertex 
V2 enters only the singlet channel, see Eq. (2.17). As we 
will see, in integrals for physical quantities containing the 
singlet channel, the main contribution comes from inte- 
gration over the momenta pi and p2 well separated from 
each other on the Fermi surface. The contribution from 
small |pi — P2I is not singular and is small. So, for di- 
mensionalities d > 1, we do not need taking special care 
about this region when performing calculations for the 
singlet channel. This argument cannot be used in d = 1 
but, in this case, the region of the momenta all close to 
each other does not give logarithms and the contribution 
coming from this region can be neglected anyway. 
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a) 



-P,(T2 



P,ai 



-P-K,ai 




P + K,a2,, 



-P-Ko 



-P-K, 



FIG. 3: a) Second order logarithmic contribution due to the 
Cooper channel, ai 7^ 02- b) The same contribution ob- 
tained as an interaction of the spin modes generated by ver- 
tices Fig. lb, and k -C Pf- This diagram coincides with the 
rcnormalization of the quadratic and quartic parts of the spin 
wave Lagrangian, see Figs. 9 b), 6a) re- written in terms of 
electronic lines. 



Finally, we notice that singling out the small and large 
angle scattering amplitudes performed here is very simi- 
lar to the one carried out for one dimensional systems^^. 
In Id the vertices Vi and V3 are indistinguishable at all. 



B. Hubbard-Stratonovich transformation 

Having reduced the interaction Cint, Eq- (2.7), to the 
interaction Smt, Eq. (2.15), we can decouple the quartic 
term by integration over an additional field 4>(t,<t' {x\ n) 
slowly varying in space. Wc will still use the notation 
(2.2), and the unit vector n labels the direction of the 
momentum on the Fermi surface. We introduce a matrix 
operator $ by defining its action on any function r]a {x) 



(2.21) 

Then we represent the 2x2 matrix in the spin space 
1 {x; n) in the form 



(j) {x, n) = iip {x, n) + cr • h {x, n) , 



(2.22) 



where Ig- is the 2x2 imit matrix in the spin space, 
(p (x, n) is a real function, h (x, n) is a three dimensional 
real vector function h = [h^, hy,hz), and cr is the vector 
of the Pauli matrices, as it was defined after Eq. (2.14b). 
The auxiliary fields are defined as bosonic: 



i(r,r; n) 



r,T-|- ^; n 



(2.23) 



Then, the partition function Z, Eq. (2.3), can be re- 
written as a functional integral over the smooth fields: 



Z=-^J z[^Ws {cl>} Wt {h} £»^Dh 
Zst = {<!>} Wt {h} DipDYi. 



(2.24) 



Here, the functional ■^|'^| is the partition function of 
the non-interacting fermions subjected to the smooth 
field 

Z[^ = j exp j dxCeff {<^}) DxDx*. (2.25) 

The Lagrangian density Cgf f for a given configuration of 
the fields (p, h reads 



(2.26) 



where the Largangian density for the fermions Lq is de- 
fined in Eq. (2.5). 

The weights for the bosonic fields Wst are of the form 

(a;, n) rfa;dn| , 

(2.27a) 



exp 



1 



'•P {x, n) 



V7V 



Wt = expl-^ I hi {x, n) 



(a;, n) dxdn 



t—x,y,z 

(2.27b) 

Here we use the notation (2.2) and the convention (1.2) 
for the integration over the direction of the momentum 
n. 

The operators 'Vs,t in Eqs. (2.27a, 2.27b) are defined 
by its action on any function a{x,n) 



t 



I 



{x,n) = ± / dri/(r-ri) 



t 



(ri,T;n) , 

(2.28) 

where Vg^t are defined in Eq. (2.18). Different signs for 

the singlet (+) and triplet (— ) channels correspond to the 
fact that the operator Vg is the positive definite and Vt 
is the negative definite, see Eq. (2.17). The choice of the 
weights (2.27) together with the fact that the fields Lp, h 
must be real is guided by the requirement that the func- 
tional integrals to be absolutely convergent. The func- 
tion / (r) is the Fourier transform of the function / (k) 
defined in Eq. (2.13): 

/(r)=ro-7o(r/ro) (2.29) 



8 



and fo{r) is the Fourier transform of fa{k). The function 
/ (r) tends to a constant r^'^'Jo (0) in the Umit |r| /ro 
and vanishes in the limit |r| oo. The role of this 
function is to regularize the theory at small distances 
leaving all interesting long distance physics intact. 

Notice that Eqs. (2.27)-(2.28) are local in time and, 
therefore, the factors Zs^t of Eq. (2.24) are not relevant 
for the determining the properties of the system. Those 
factors will be usually suppressed in the subsequent for- 
mulas. 

Thus, we have decoupled the interaction term Sint , Eq. 
(2.15), in the action S, Eq. (2.4), with the Hubbard- 
Stratonovich transformation, Eqs. (2.24)-(2.27b). As 
the field (j){x,n) varies slowly in space, we can apply 
quasiclassical description for the electron Green functions 
in such fields. 



C. Quasiclassical Green functions 

Let us introduce Green functions Ga,a' (a;, |(^|^ 

corresponding to the Lagrangian density Cgff > Eq. 
(2.26), as 

G„,y{x,x'\{^})=Z-'{$} (2.30) 
X Jxa {x)x:' {x')exp (^-J dxCeff [$]) DxDx*. 

In what follows we will suppress the argument | |(^| 
whenever its presence is self-evident. 

The Green function G ^a;, a;'| li^l) : is the matrix 2x2 
in the free space, is the functional of real fields ip, h, and 
it satisfies the equations 

-^-Hor^G + ^G = 6{x-x')la, (2.31a) 

^ - ^or') G + G$ = 5{x-x') K, (2.31b) 

where Hor, Eq. (2.6), acts on r. The action of the opera- 
tor $ on the Green function is determined by Eq. (2.21) 
as 



$G 



G^ 



^<P<7,a" ( j^2^ ''^' |~|^ G^//<^/(ri,r;a;')| ; 

X |G<^(7"(a;;ri,r')?!'<7",<7' ( j^ '^^ ,t'; |~| ^ | • 

(2.31c) 

The derivation of the equations for the quasiclassical 
Green functions can be carried out in the same way as 



in Rcf. 48. Subtracting Eq. (2.31a) from Eq. (2.31b) and 
making a Wigner transformation 

G (^; ^0 = / -^^'""G {r, r'; p) (2.32) 



of the result, we obtain 
d d «pVj 







+ 



dr dr' 



m 



G(T,T';r,p) 



(2.33) 



G (r, r'; r, p) ^ (r, r'; n) - 4> (r, r; n) G (r, r'; r, p) 



Equation (2.33) is justified provided the dependence 

of the field (j){x,n) and, hence, of G(t, T';r;p) on the 
coordinate r is slow on the scale of the order of Fermi 
wavelength. This is guarded by the cutoff scale ro in 
Eq. (2.13). In principle, one could derive Eq. (2.33) more 
accurately, which would produce additional terms con- 
taining phase space derivatives of the functions (j) and G 
in the second line. However, the additional derivatives 
would suppress the infrared singularities we are inter- 
ested in, and that is why we neglected them. At the 
same time, no higher derivatives arise in the first bracket 
in Eq. (2.33) and this term is exact for the quadratic 
spectrum of the fermions, Eq. (2.6). 

The next step is to reduce the Green function G(p, r) 
to a function involving the degrees of freedom describing 
the motion of the system along the Fermi surface. To 
accomplish this task we linearize the spectrum by putting 



P 

m 



(2.34) 



in Eq. (2.33), where n is the unit vector. The justifica- 
tion of such approximation is that the fluctuating fields 
can mix the electron states only in the vicinity of the 
Fermi surface whereas the states deep in the Fermi sea 
remain intact. After approximation (2.34) the operators 
in Eq. (2.33) do not depend on the variable 



£, = p^/2m - ep, 



(2.35) 



that describes the evolution perpendicular to the Fermi 
level and the latter can be integrated over. Then, one 
obtains*^'^^ 

^ + ^-'^^^n-pVrj5(T,T';r,n) (2.36) 
+ (T,r';r,n)0(r,T';n) - (^(r,T;n)5(r,r';r,n) , 



where 
g(T,T']r,n |(^|) 

G 



T,T';r; {PF + — 
vf 



n|{<^} 



(2.37) 



The function g must obey the antiperiodicity condi- 
tions 



g(r,T';r,n) = -ff(T+l/r,r';r,n) 
= -g{r,T' + 1/T;r,n) 



(2.38) 
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that follow from Eqs. (2.8), (2.30), (2.32), and (2.37). 
Clearly, this condition is consistent with Eq. (2.36) for 
the periodic fluctuating fields see Eq. (2.23). 

Equation (2.36) is linear and therefore is not sufficient 
to find g (t, r'; r, n) unambiguously. In order to define the 
problem completely, one has to complement Eq. (2.38) 
with a certain constraint. To derive this constraint, we 
introduce a new function 

s(T,T';r,n|{0}) 

/■l/T 

= / « (t,t"; r,n 5 (r", r'; r, n [4>]) . 

J 

(2.39) 

Using the definition (2.39), we obtain from Eq. (2.36) 







+ 



^ + ^-«^i^n-pVr)-B(T,T';r,n) (2.40) 



B (r, r'; r, n) ^ (r, r'; n) - 4> (r, r; n) 13 (r, r'; r, n) 



In the absence of the fluctuating flcld, (/) = 0, the Green 
function can be easily found from Eqs. (2.31a), (2.32) and 
(2.37) 



»(.,/;r.n|0) = ii^rj:.M.-.,y_i 



Re 



sin7rT(r-T' + iO) 



(2.41) 



Substituting Eq. (2.41) into deflnition Eq. (2.39) and per- 
forming the integration, we find 



S (t,t'; r,n|o) = (r - r') ■ 



(2.42) 



Substitution of Eq. (2.42) into Eq. (2.40) shows that 
Eq. (2.43) remains valid for the arbitrary field (j), 



B(T,T';r,n|{0}) =1,5(t-t'). 



(2.43) 



i.e. no perturbation by the Hubbard-Stratonovich fields 
can violate the condition (2.42). 

Equations (2.39) and (2.43) complement Eq. (2.36) 
and these equations are sufficient to find the function g. 
Eq. (2.36) is much simpler than the original Schrodinger 
equation (2.31a) as it operates with the smooth quanti- 
ties and involves only the first derivatives. The further 
program is to solve Eqs. (2.36), (2.39) and (2.43) for ar- 
bitrary configurations of the fields 4>- After that, in or- 
der to calculate physical quantities, one should perform a 
proper averaging over fields (j) with the weights defined in 
Eqs. (2.27). All this is still not a simple task and in the 
next Sections we will express the solution of these equa- 
tions in terms of a functional integral over supervectors, 
in order to obtain the local theory in terms of only the 
bosonic variables describing the collective excitations. 



III. CHARGE AND SPIN COLLECTIVE 
VARIABLES. PARTITION FUNCTION 

A. Further simplification of the quasiclassical 
equations 

Solutions of Eqs. (2.36), (2.39), and (2.43) describe col- 
lective excitations and our task is to find them at least 
symbolically in order to facilitate calculation of the par- 
tition function Z{(j)}, see Eq. (2.25) and the averaging 
over the auxiliary field (j). 

At first glance, we coiild simply follow the scheme de- 
veloped in Ref. 35 writing the solution of these equations 
in terms of a functional integral over constrained super- 
matrices. However, in the present situation this scheme 
is not convenient due to the dependence the Hubbard- 
Stratonovich field (j) on T. 

Instead, we look for the solution of Eqs. (2.36), (2.43), 
and (2.39) in a form 

g (r,r'; r,n| j^j) 

= t (r; r, n| {,^}) g (r, r' |o) f-' (r'; r, n| {^}) 

(3.1) 

where the Green fimction for the free electrons, 
t,t' 0^ is defined in Eq. (2.41). The 2x2 matrix 

in the spin space, T satisfies the condition 

t (r; r, n| {<^}) = t (^r + i; r, n| {<^}) (3.2) 

so that the antiperiodicity of the Green function (2.38) 

is preserved. In the remainder of this subsection, we will 

suppress the argument | j^j whenever it is self-evident. 

The representation of the Green function in the form 
of Eq. (3.1) is nothing but the matrix form of the eikonal 
approximation, which can also be viewed as a generaliza- 
tion of the Schwinger Ansatz^*^. It easy to check that the 
Green function, Eq. (3.1) is consistent with Eqs. (2.43), 
and (2.39), and what remains to be done is to find the 
proper matrix T, such that Eq. (2.36) is satisfied. 

Substituting Eq. (3.1) into Eq. (2.36) we obtain 



5(r,r'|o) [KiT,v,n)-K{T',T,n) 



0, 



(3.3) 



where the 2x2 matrix in the spin space K is given by 

K (x, n) = {x, n) {dr - ivpnVr) f (x, n) (3.4) 

— {x, n) (j) {x, n) T {x, n) , 

and we use the short-hand notation (2.2). 

Equation (3.3) must be fulfilled for any t and t'. This 
is possible only for drK^XjU) = 0. Using Eq. (3.4), we 
obtain 



{—dr + zuFnVr) T {x, n) 

= T {x, n) A (r, n) - 4){x, n)T {x, n) 



(3.5a) 
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where A (r, n) is an arbitrary time independent matrix. 

We can transform Eq. (3.5a) to a more convenient form 
writing the corresponding equation for (x, n) 

{dr - ivpnVr) {x, n) (3.5b) 

= A (r, n) (x, n) - f - 1 (x, n) 4){x, n) 

We differentiate Eq. (3.5a) with respect to r and 
post-multiply it by {x,n). Then we pre- multiply 

Eq. (3.5b) by (x, n) and subtract thus obtained equa- 
tions from each other. As the result, we find 



where 



{—dr + ivpu.'Vr) M (x, n) 

+ (j) {x, n) , M„ (x) = -dr${x, n) 

df{x,n). 



(3.6) 



M (x, n 



dr 



-r-'{x,n) 



(3.7) 



and the symbol [...,...] stands for the commutator. 

Using the representation (2.22) for the matrix (j) (x, n). 
we look for the matrix M (x, n) in the form 



M (x, n) = ip (x, n) + S (x, n) • cr 



(3.8) 



where p (x, n) is a scalar real field and Sn (x, n) is a real 
three dimensional vector field. As follows from Eqs. (3.8) 
and (3.2), those fields are periodic 



P (t, r, n) 



p ( T+ -,r,n 



S (r,r,n) = S ( T + -,r,n 



(3.9) 



Substituting Eq. (3.8) into Eq. (3.6) we obtain two 
independent equations for p(x,n) and S (x,n): 



-I- W^nVr ) Pn {x) 



dr 



dip (x, n) 



(3.10a) 



+ iu^nVR Sn (x) 



+ 2i[hn (X) XSn (X)] =- 



(x) 
dT 



(3.10b) 



It is easy to see that Eqs. (3.10) are consistent with the 
periodicity requirements (3.9) and (2.22). 

Equations (3.10) are the final quasiclassical equations 
that will be used for further calculations. We emphasize 
that Eqs. (3.10) are obtained from Eqs. (2.36). (2.39). 
(2.43) without making any further approximation. The 
field p (x, n) corresponds to the density fluctuation in the 
phase space, whereas the field S (x, n) describes the spin 
fiuctuations. 

Equations (3.10a) and (3.10b) determining these fluc- 
tuations due to the Hubbard-Stratonovich fields arc re- 
markably different from each other. Equation (3.10a) for 



the density is rather simple, and can be solved immedi- 
ately by the Fourier transform. This is what one obtains 
using the high dimensional bosonization of Refs. 6,7,39- 
46 from an eikonal equation. Of course, we could take 
into account gradients of the field (x, n) and this would 
lead to additional terms in the L.H.S. of Eq. (3.10a). 
However, this does not lead to new physical effects. 

In contrast, Eq. (3.10b) is not readily solvable due to 
the presence of h (.r, n) in the left-hand-side (L.H.S) of 
this equation. Actually, the L.H.S of Eq. (3.10b) is just 
the equation of motion of a classical spin-density in the 
external magnetic field h. We will see that the presence 
of this form will result in non-trivial effects that will be 
considered later. To the best of our knowledge, this dif- 
ference between the charge and spin excitations in > 1 
has not been emphasized in literature. 



B. Partition function 

Having found the semiclassical representation for the 
Green functions, we are prepared to express the partition 

function j*^! from Eq. (2.25) in terms of the collective 

variables p (x, n) and S(x, n). Integrating over x, in 
Eq. (2.25) and using Eqs. (2.26), (2.5) for the Lagrangian 

density >Ce// l*^!; write l*^} i^i the form 



InZ |(^| =Tr j\n (^-drK - HqK 



(3.11) 



where the Hamiltonian Hq is introduced in Eq. (2.6), 

operator $ is defined by Eq. (2.21), and Tr includes the 
trace in the spin space as well as the integration over r, r. 

Equation (3.11) can be rewritten using the standard 
trick of integration over coupling the constant as 



InZj^l -lnZ{0} 



duduTr / In -drK - Hota + 



) 



J du"^ J dx l>G {x,x), 



(3.12) 

where the Green function is obtained from 

that of Eq. (2.31a) by the rescaling of the Hubbard- 
Stratonovich fields: $ — > and the action of the op- 
erator $ is defined by Eq. (2.31c). The term lnZ{0} 
describes the thermodynamics of the non-interacting 
fermions, and we will suppress this term in all the subse- 
quent formulae. 

Using Eqs. (2.31c) and (2.32), we obtain from 
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Eq. (3.12) 
In Z 

= J du J dx J -^^Tra^) {x; n) G {t, t; r, p| {m^}) , 

(3.13) 

where Tr^- denotes the trace in the spin space, and the 
short hand notation (2.2) is used. We represent the inte- 
gration over the momentum as 

where i/(^) is the density of states (DoS) per one spin ori- 
entation, ^ = corresponds to the Fermi level and we use 
the convention (1.2) for the integration over the direction 
over the momentum on Fermi surface, n. Neglecting the 
energy dependence of DoS, we obtain from Eq. (3.13) 

lnz|(^| 

1 

= —i-KV j du j j dnTra4> {x; n) g ^r, r; r, n| |"<^|) i 

(3.14) 

where v = ^{^ = 0) is the DoS on the Fermi level per 
one spin orientation, Green function g ^r, r'; r, p| 
satisfies the constraints (2.39) and (2.43) and satisfies 
Eq. (2.36) with the rescaling (f) u(j). According to 
Eq. (2.41) g is a singular fimction at coinciding time, 
so that the equal time value should be understood as 

g{T, t) = \ lim [g{T, t + 5) + g{T, t - S)] . (3.15) 

Next, we substitute Eq. (3.1) into Eq. (3.14). Using 
Eq. (2.41) and the rule (3.15), we find 

g (r,r; r,n| {w^}) 

t^r; r,n ^^u$^^ f'^ (t + 5;r,n [u4>]) 

■ f (t; r, n| [u^Y) " ^'^ ^' "| {^^}) 



Here 



— lim — 

2 5^0 6 



dT 



r; r. n 



Using the definition, Eq. (3.7), the representation (3.8), 
and equations of motion (3.10), we obtain finally from 
Eq. (3.14) 

=ZoZpMZ,{h}. (3.16a) 



Zp = exp 



Zs = exp 



2z/ y du J p {x, n; u) (p {x, n) dxdn 



(3.16b) 



-2^ J du J S {x,n;u)h{x,n) dxdn 



(3.16c) 

The functions p {x, n; u) and S {x, n; u) should be found 
from the equations [cf. Eqs. (3.10)] 



{—dr + WftoVt) P {x, n; u) = —udr(p {x, n) , (3.17) 

LuS {x,n;u) = —udrii{x,n) . (3.18) 
In Eq. (3.18), the operator Lu equals 

Ly, = {-dr + ivpnVr) + 2iuh{x, n) (3.19) 
where the matrix h has the following form 



—hz{x,n) hy{x,n) 
h{x,n) = I hz{x,n) —hxix,n) 

-hy{x,n) hx{x,n) 



, (3.20) 



and hx,hy, and hz are the components of the real vec- 
tor h (ha = [h X a] for any vector a). We will call this 
space of three dimensional vectors the "spin space" as 
the 2x2 spin space for the original electron will be no 
longer needed in further considerations. The functions 
S{x, n; u) and p{x, n; u) satisfy the periodicity condition 
(3.9). 

The operator L„, Eq. (3.19), is antisymmetric 



L,, — — £/,, 



(3.21) 



where the transposition "T" includes both the changing 
of the sign of the derivatives and the transposition of the 
spin indices. However, this operator is neither Hermitian 
nor anti-Hcrmitian. The importance of this subtlety will 
be underlined in the next section. 

Thus, in order to calculate the partition func- 
tion, Eq. (2.3), for the system of interacting fermions 
in the quasiclassical approximation, one should solve 
Eqs. (3.17), (3.18) and substitute their solutions into Eqs. 
(3.16a)-(3.16c). Then, one should use Eq. (2.24) and av- 
erage over the fields (p and h with the weight given by 
Eq. (2.27). 

Before proceeding, we notice that there is a well-known 
flaw in the quasiclassical approximation (3.14) to the ex- 
act Eq. (3.12) (this flaw is usually referred to as an ul- 
traviolet anomaly). In Eq. (3.12), the two times in the 
Green function are put equal to each other before the 
integration over the momentum is performed, whereas 
Eqs. (3.14)-(3.15) imply the opposite order of the limits. 
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Those operations do not commute as they treat contri- 
butions from the region far from the Fermi surface differ- 
ently: in the quasiclassical approximation the informa- 
tion that the electron states are limited from below at 
^ > —€f is lost. Lost contributions, however, are coming 
from the transitions with the large energy and therefore 
are perfectly analytic functions of fields ip, h, and their 
gradients. As the result, Eqs. (3.16a) is modified as 

Z {0} = ZoZ, M {^} Z, {h} ZT {h} 

InZ^^ = c^a;rfnrfnXa;, n)^^ (rm') 4>{x, n') + • • • 

InZ™ = y rfa;rfnrfn'h(a;, n)4 (^rm'^ h(a;, n') . . . , 

(3.22) 

where . . . stand for the terms containing higher gradients 
of the field or the higher powers of the field. All such 
terms however will be small as l/ci? and that is why 
keeping them would be the overstepping of the accuracy 
of the quasiclassical equations (2.36). 

Functions dp^g (9) depend on the details of the ul- 
traviolet cut-off [for the weakly interacting gas dp = 
da]- One property, however, remains intact - the 
response of the system on the fields independent on 

the coordinate but arbitrary periodic function of time 

\ It \ It 

(y9(T), h(T); /p ip{T)dT = /p h(T) = should vanish, 

because the total charge and the total spin commute with 

the Hamiltonian: 

Zp {^{t)} Z7 {^(t)} = Z, {h(r)} Z7 {h(r)} = 1. 

(3.23) 

For such fields, Eqs. (3.17) - (3.18) are trivially solved 
p(r, u) = —u(j){T), S(r, u) = — uh(T) and we obtain from 
Eqs. (3.23), (3.16b), (3.16c) and (3.22) 



dp ( nn' 



dw ■ 



da I nn' I dn = 1. 



(3.24) 



All the other properties of dp^s are model dependent and 
can be established by direct perturbative calculation for 
stationary fields for which the semiclassical contributions 
(3.16b)- (3.16c) vanish. 

However, it would be a redundant exercise, as Eq. 3.22 
have the same form as the weights (2.27), and the role of 
terms (3.22) is just a renormalization of the constants in 
those weights. The contribution of the interaction terms 
with the high-momentum transfer not included into the 
Hubbard-Stratonovich transformation leads to the simi- 
lar effects. It means, that the form of the weights for the 
fields h, should be established not from the first princi- 
ples but from the requirement that the quadratic part of 
the theory should reproduce the bosonic modes obtained 
from the kinetic equation in the Landau theory of Fermi 
liquid. It leads to the replacement of Eqs. (2.27) with 



W.^exp{-^/„.,„,[f;V](x,n),t.<in} 



Wt = exp <^ --^ ^ \ h, (x, n) 



(x, n) dxdn 



(3.25b) 

where the operators Ts,t are defined by its action on any 
function a {x, n) as 



2r, = / 

2ft = -/ 



1 -l-F^ 



(3.25c) 
(3.25d) 



Here operators / and F are defined by their action on an 
arbitrary function 6(t, r; n) as 



fb 



(r,T;n)= j dri/ (r - ri) 6(ri, r; n) 



(3.25e) 



(r, r; n) = / dnzF"''^ (nna) 6(r, r; na). 



and the convention (1.2) is used. The cutoff function 
/ is defined by Eq. (2.29), while the functions ¥P''^ are 
the Fermi liquid functions describing the interaction be- 
tween two quasiparticles in the singlet or triplet states. 
We will see in the next subsection that the choice (3.25) 
indeed reproduces the correct propagators for the collec- 
tive modes in the Fermi liquid theory. In what follows 
we assume the operators Ts,t to be positive definite, and 
the system far from Pomeranchuk instabilities. 

With the help of the quasiclassical consideration we 
can recast Eq. (2.24) to the form 



n = -T\nZ = Oo + fip + 



(3.26a) 



(3.25a) 



where fio — —TIuZq describes the leading contribution 
of the quasiparticles. The leading singular corrections are 
associated with the collective modes and they are given 
by 

exp(^-^) = J D^WsW}Zp{v>} (3.26b) 

exp(|-^) = J DhWt{h}Z,{h}, (3.26c) 

where the functionals g are given by Eqs. (3.16b)- 
(3.16c). In writing the expression for the partition fimc- 
tion we ignored the terms, [e.g. Zgt in Eq. (2.24)], which 
do not lead to change of any observable quantities. 

The results of the present Section show that study of 
the system of the interacting fermions can be reduced 
to investigation of a system of bosonic charge and spin 
excitations. Therefore the word "bosonization" is most 
suitable for our approach. We see that the method should 
work in any dimension. At the same time, it is more gen- 
eral than the scheme of the high dimensional bosoniza- 
tion of Refs. 6,7,39-46 because we can consider the 
spin excitations that are much less trivial than the charge 
ones. 
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C. Thermodynamics of free modes. 

Before we start formulating the proper field theory de- 
scription for calculating the partition functions (3.26b) 
- (3.26c), it is instructive to try to determine it by the 
brute force analysis of Eqs. (3.16b)-(3.16c). 

For the charge mode we immediately solve Eq. (3.17) 
by the Fourier transform: 



p{uj„,k;n) = u- 



lUlr, 



ip{uJn,k;n). (3.27) 



where u>n = 27rTn is the bosonic Matsubara frequency, 
Eq. (3.9), and 



93(w„,k; n) = 9?*(-w„, -k; n) 



(3.28) 



because the field (p{x; n) is real. Substituting Eq. (3.27) 
into Eq. (3.16b), we find 



Zp = exp 



d'^kdn iujn |</?(w„.k; n) 



lU)n 



Uirkn 



(3.29) 



Both the partition function, Eq. (3.29), and the weight, 
Eq. ((3.25a)) are the Gaussian functionals and, therefore, 
the fimctional integration in Eq. (3.26b) can be readily 
performed with the result 



T rd^ 



l + /(k)F^- 



^ + VF^n 

(3.30) 

Here, factor of 1/2 originates due to the constraint (3.28), 
and the action of the interaction function is defined 
by Eq. (3.25e). The function / (k) , Eq. (2.13), cuts 
momenta k exceeding r^^ ■ However, its presence in Eq. 
(3.30) is important only for calculation of corrections to 
the coefficient in the linear term in the specific heat. Non- 
trivial contributions to the specific heat C come from the 
momenta k ~ T/vp r^^ and do not depend on the 
function /. The explicit formulas for the specific heat 
given by Eq. (3.30) will be derived in Sec. VII, how- 
ever, here we will give the equivalent representation of 
Eq. (3.30) more convenient for the comparison with fu- 
ture material: 



d'^kdn . 
7^ 



In 



l + /(k)7^ 



tOJ„ 



r 



iij^kn 
(3.30') 



Equation (3.30) has a very simple form and de- 
scribes the thermodynamic potential of the collective 

non-interacting charge mode. In the conventional dia- 
grammatic language, Eq. (3.30) corresponds to the con- 
tribution of ring diagrams. The advantage of the deriva- 
tion here is the explicit demonstration that Eq. (3.30) 
completely solves the problem of the singular corrections 



in the charge channel [which is the only one present for 
the spinlcss electrons]. No further terms are present for 
the linearized spectrum and all of the other corrections 
have additional smallness of the order of T/ep in com- 
parison with Eq. (3.30). This means that no further con- 
sideration of the singlet channel is necessary. 

Let us turn to the triplet channel. One can see that 
due to the presence of the Hubbard-Stratonovich field 
h in the operator (3.19), one can solve Eq. (3.18) only 
approximately. In particular for |h| —^ 0, one finds [cf. 
Eq. (3.27)]: 



S {tOn, k; n) = u- 



IU!„ 



iuin — Virkn 



h{ujn,k;n) + (3.31) 



where . . . stand for the functionals of the second and 
higher orders in field h, and 

h(a;„,k;n) = h*(-a;„,-k;n). (3.32) 

Substitution of Eq. (3.31) into Eq. (3.16c) yields 



exp 



vT 



lUJri 



d'^'kdn 



j=x,y,z 



\hj{u)n.k;n)\ +... 



(3.33) 



where . . . denote the functional of the third and higher 

orders in h. It is important to emphasize that such non- 
linear terms do not have any additional smallness inT/ep 
in contrast to the singlet channel formula (3.29). 

If we ignore those non-linear terms we end up with the 
Gaussian functional integral in Eq. (3.26c) and we obtain 
analogously to Eq. (3.30): 



l.(°) = ^V^ln 



l + /(k)F- 



vphn 



-iu)n + Vp^in 

(3.34) 

and the action of the interaction function F'^ is defined 
by Eq. (3.25e). The additional factor of 3 in comparison 
with Eq. (3.30) stands for the three independent compo- 
nents of the spin density. The equivalent representation 
for Eq. (3.34) is [cf. Eq. (3.30')] 



0„ 



3r 

~2 



E 



d'^Mn . 
7^ 



In 



l + /(k)7 



lUln 



Vp\i\\ 



ILOr, 



^-FF<^' 



Vplin 
(3.34') 



Let us mention for the future comparison with previous 
works, that for the weakly interacting systems, the ker- 
nels 7 are the linear function of the scattering amplitudes 
(2.17) 



lp = '^Vs; ^ = -vVf 



(3.35) 
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However, due to the presence of the ignored non-smaU 
terms in Eq. (3.33), formulas (3. 34), (3. 34') are by no 
means exact or correct low temperature asymptotic ex- 
pression. The next section present an efficient calcula- 
tional scheme to deal with this non-linearity. 



IV. SUPERSYMMETRY APPROACH FOR THE 
BOSONIC EXCITATIONS. 

As we have already explained, the exact derivation of 
the functional {h} determining the thermodynamics 
of the triplet mode Eq. (3.26c) is not possible. Moreover, 
obtaining the non-linear terms by the further expansion 
of the solution in powers of h is not a correct way to pro- 
ceed because the resulting theory in terms of h will be not 
only non-linear but also non-local which would obscure 
such important features of the theory as its renormaliz- 
ability. 

An analogous problem exists in the theory of disor- 
dered metals but in many cases it can be overcome using 
the supersymmetry method^^'^^. The main idea of the 
method is to express the solution of a linear equation with 
a disorder in terms of a functional integral over auxiliary 
supervectors containing both conventional complex num- 
bers and anticommuting Grassmann variables. Then, one 
is able to average over the disorder and reduce the dis- 
ordered system to a regular model without any disorder 
but with a local effective interaction. 

We will borrow these ideas. The role of the disorder 
here will be played by the field h itself, so its imaginary 
time dependence will lead to certain modifications. 

We will discuss the number and the properties of the 
necessary fields in Subsection IV A in somewhat simpler 
form and will write down the full-fledged effective La- 
grangian in Subsection IV B. The final form of the field 
theory obtained after the averaging over field h is given 
in Subsection IV C and it is generalized further in Sub- 
section IV D. 



A. The number of auxiliary fields and 
Hermitization. 



Solution of Eq. (3.18) can still be written in a symbolic 
form as 



S {x, n; u) = —uL^ ^drh {x, n) 



(4.1) 



where the local operator L„ is given by Eq. (3.19). Sub- 
stitution of Eq. (4.1) into Eq. (3.16c) yields 

Zg = exp 2v J udu J dxdnh{x,n)L~^drh{x,n) 

(4.2) 

and we use the notation (2.2) throughout this section. 

The argument of the exponent is non-local and our goal 
is to get rid of such a non-locality. The standard route 



to proceed would be to re-write Eq. (4.2) as a functional 

integral 

{h} = J DS^'DSDxDx' exp - '^'^ J du j dxdn 
X {S*L„S + x*LuX + uSdrh + S*h} 



(4.3) 



Here 




(4.4) 



The spin space s was introduced after Eq. (3.20). The 
fields S*, S are the usual complex vector fields and x, X* 
are anticommuting Grassmann fields needed to cancel out 
the operator determinant. All the fields are functions of 
X, n, u and satisfy the periodicity conditions ^\ 



1 



(4.5) 



S(a;,n,u) = S ( r,r+ -,n,u 



S* {x, n, u) = S* ( r, r + -, n, u 



Xix,n,u) = x{r,T+ -,n,u 



X* {x, n,u)=x* {r,T+-,n,u 



Because the Grassmann fields are periodic rather than 
anti-periodic we will call them pscudofcrmions. The ar- 
gument of the exponent (4.3) would be, then, the local 
functional linear in h, so the integration in Eq. (3.26c) 
could be easily performed yielding the local expression in 
terms of powers of S and % only. 

However, Eq. (4.3) is rather deceptive. Indeed, the 
possibility to write such a functional integral is based 
on the assumption that the integration over the bosonic 
fields is defined for an arbitrary configuration of the field 
h, in particular the directions of such integrations cannot 
depend on the field h at all. Therefore, Eq. (4.3) would be 
correct only if the operator !,„ of Eq. (3.19) were positive 
definite which is not the case. Moreover, as we have 



Note that although the theory is still Gaussian we had to intro- 
duce one more coordinate u, such that the fields depend on u for 
the representation of the needed determinants. The other possi- 
ble way would be to use representation (3.7)-(3.8) for S and write 
WZNW action for matrices T, see e.g. Ref. 51. We chose not 
to pursue this line because of the difficulty of the identification 
of the manifold of the integration over the matrices T such that 
the "saddle point" matrix T solving Eq. (3.5b) would belong to 
the integration manifold. 
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mentioned after Eq. (3.21), L-u is not even Hermitian and 
we do not have a priori the knowledge about the signs of 
real and imaginary parts of eigenvalues. 

Those complications make expression (4.3) mathemat- 
ically meaningless. Fortunately, the proper procedure 
for non- hermitian operators containing first-order deriva- 
tives has also been worked out prcviously^^ and we will 
use this method for further calculations. 

Let us double the number of the bosonic and pseud- 
ofermionic fields as 



X 



H 



X' 



H 



(4.6) 



st = ([si]*,[sT)^; x^ = {[xr,[x'r)H, 



where each element has the structure of Eq. (4.4). We 
will call the additional space of the two- component vec- 
tors "hermitized" space and use the subscript {■ ■ ■)h for 
the writing explicitly the structure in this space. We 
define the new operator M„ acting in this doubled space 



(4.7) 



where each element of this matrix is a matrix in the spin 
space, sec Eq. (3.19). By construction, the operator, Eq. 
((4.7)), is Hermitian, M = M\ and thus the functional 



j dx (s^M„s) 



is real. Therefore, the identity 



1 = j DS^DSDxDx^e^v[-Sf] 
-2iy j du j dxdn 

X {st (m„ + ts)s + (m„ + id) x} 



Ceff 



holds for an arbitrary configuration of the field h as the 
integral over the bosonic fields is always convergent. Here 
(5 is a positive real number and the limit S +0 is to 
be taken at the end of the calculation. The fields satisfy 
the periodicity conditions, Eq. ((4.5)). Finally, using 
Eq. (4.8) and obvious formula 



1 



Lu 
1 



1 

TT' 
1 



ft ' 



1 

J-'u 

1 1 

J-'u J-'u ' H 



we obtain instead of Eq. (4.3) 

{h} = j DSDSDxDxexp [S^' 

v^2i J du J dxdn 



exp 



{u (S^ + S^) drh + ( 



>} 



(4.9) 



Equations (4.8) (4.9) is the final result of this sub- 
section. We have succeeded in rewriting the original 
non-local expression written in terms of the Hubbard- 
Stratonovich field h into the theory local in terms of the 
new fields S, %. As this action is a linear functional in h, 
we will be able to integrate it out and obtain the action 
in terms of those new fields S,x only. Before doing so, 
we will recast Eqs. (4.8)-(4.9) in a more compact form. 



B. Supervectors and the effective Lagrangian. 

We introduce the superspace (graded space) as the 
space of vectors having the same number of complex and 

Grassmann components''^. In particular, the field intro- 
duced in Eq. (4.6) can be compactified as one supervec- 
tor^" 



(4.10) 



where subscript g stands for the graded space, and each 
element has the structure of Eq. (4.6), so the fields are 
defined in a linear space obtained as a direct product 
of Hermitized (H) spin (s) and superspace (g). In other 
words, notation ip means the 12 -component supervector 
defined in a linear space s H g. 

Using notation (4.10) the action from Eq. (4.8) can be 
re-written in a short form 



(4.8) "^h*^ = -2«z^ J du J dxdn ip^ [Mu 



(4.11) 

where Ig is the 2x2 unit matrix acting in the superspace. 

We will see shortly that the most interesting contri- 
bution will come from the scattering terms where the 
direction of the spin changes its direction to opposite. 
Anticipating this fact, we will join ^'(n) and y>(— n) in 
one vector of the larger dimensionality 

'^(-) = (vK-n))j -^'("^ = (vnn);V't (-n))^, 

(4.12) 



^" Using the same notation for the supervector here and for the 
Hubbard-Stratonovich field in Sec. II should not lead to a con- 
fusion as the latter will not appear in any further consideration. 
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where each clement has the structure of Eq. (4.10). and 
keep integration over n in all of the subsequent formulas 
over the d-dimensional hemisphere, say Ux > 0. The final 
answers definitely will not depend on particular choice of 
the hemisphere. From now on the integration over the 
momentum direction will mean 



dn - ■ ■ 



dn. . .: 



nx>0 



Idn = -. 
2 



(4.13) 



We will call two-dimensional space defined in Eq. (4.11) 
the "left-right" space, using the analogy with the one- 
dimensional systems and will denote it by subscript n. 

The final step in the definition of the supervector is 
once again performed for the calculational convenience, 
and it is equivalent to the introduction of the Gorkov- 
Nambu spinors in the theory of the superconductivity^^. 
We increase the size of the supervector as 



1 

71 



m*;<P*).„ (4.14) 



where each component has the structure of Eq. (4.12). 
We will call the corresponding space "electron-hole 
space" and denote it by the subscript eh when written 
explicitly. The benefit of this doubling of the size of the 
supervector will become apparent when the perturbation 
theory for the resulting model is developed in Section 
V. To avoid misunderstanding, we emphasize that the 
"electron- hole" introduced here has nothing to do with 
the electrons and holes in the original system and intro- 
duced here to label the fields of the spin excitations only. 

Using the standard convention for the complex conju- 
gate of the Grassmann variables. 



[x*]* = -X, 



(4.15) 



and definitions (4.12) - (4.14) we rewrite Eq. (4.11) as 

= -2iv J du J dxdn 'ip'' (Xiu ^tg + iS^rp ■ 

° (4.16) 
Here is the matrix in the H ®n® eh space and it 
has the structure 



Mu = 







eh 



(4.17) 



M,(n) 



M„(n) 
M„(-n) 



with matrix M„ given by Eq. (4.7). 

To make the notation consistent with the previous 
work^^'^'^'^^, we introduce the conjugated supervector as 



■0 = 1/'^ A; A = lg(g)l„(g)ls(g)le/,( 



1 
-1 



. (4.18) 



Using the explicit structure of the supervectors, Eq. 
((4.14)), and the convention (4.15) one can verify that 
the conjugated supervector ij} is related to ■0 as 



C =ts®tn' 



Co 

-Co 



Cl 



-1 

1 



C2 



eh 



H 

1 

1 



^ / Cl 
Co={ Q gj^ (4.19) 



eh 



Accordingly, the conjugation of supermatrices is intro- 
duced as 



A = CA' C 



TriT 



(4.20) 



for an arbitrary supermatrix A acting in s®g®H®n®eh 
space. For the two supervectors t/'j 2 of the structure 
(4.14) one finds 



(4.21) 



Substituting definition (4.18) into Eq. (4.16) and find- 
ing explicit form of from Eqs. (4.17), (4.7) and 
(3.19)- (3.21), we obtain 

Sf = -2iu j {X) Chip {X) dX, (4.22) 

where we use the short hand notation 
X = (r, r, n, u) ; 

J dX... = J dr J^^ dr J dnj^ , 



du.. . 



(4.23) 



and the convention (4.13) for the angular integration. 
The Lagrangian in Eq. (4.22) is given by 



Ch = Co — 2iuTs6jC.h — iSA, 



(4.24) 



where the matrix in the spin space h is defined in 
Eq. (3.20) and the free propagation Lagrangian has the 

form 



Co = ~ivF (nV)f3S3 - drAi, 
cl = -ivF (n V) fgSa + ^^Ai, 



(4.25) 



The rotation of the spin excitation by the Hubbard- 
Stratonovich field [cf. Eq. (3.20)] is described by 

/ -M^{x,n) E.y{x,n) ' 
dAh{x,n)=i H^(a;,n) -EI^(a;,n) 

\-Hy(a;,n) m^{x,n) 

I^(x,n) = lgOl^fOle„<^ ''''^^^'"^ ° 



h^{x,—n) 



H 



(4.26) 
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The supermatrices in Eqs. (4.25) and (4.26) are intro- 
duced as 



fs = 1^ O 1„ O Ijf O 1„ I 



1 

-1 



Eg = 1, O 1„ (g) Iff . 



1 
-1 



h 



(4.27) 



1 



Ai = 1^ Ig ( ^ Q 1 1„ O teh, 



H 



The action is supersymmetric, i.e. invariant with 
respect to all possible homogeneous rotations in g-space. 

To complete the derivation, we have to express the ex- 
ponent in Eq. (4.9) in terms of the supervector xp. Those 
are only terms that break the supersymmetry and thus 
lead to finite contributions to physical quantities. Using 
the definitions (4.14), (4.12), (4.6), (4.4), we find 













j du j dxdn ^ 


Si' 






> 













2 j dx(^^{X)¥\X) 



a \Vh V-J eh V'T 



/i-y(n) 
/i-y(— n) 



(4.28a) 



and, analogously. 



duj dxdn (s^ + S^) h 



1 

-1 



H 



eh 



udrhj{n) 
udrh-y{—n)^ 

(4.28b) 



where Eq. (4.23) is used, and -y = x,y,z labels the com- 
ponents in the spin space. 

Equations (4.28) and (4.22) enable us to obtain a rep- 
resentation of formulas (4.8)-(4.9) in supersymmetric no- 
tations 



Z,{h} = J £>Vexp[-5fT 

2i/V2i j '^{X)F{X)dX 



X exp 



(4.29) 



where .S^^^ is given by Eq. (4.22) and 

F{X) = F^{X) + F^{X). 



(4.30) 



The superfields in Eq. (4.29) satisfy the periodic bound- 
ary conditions 



ip (t, r, n; u) = V (r + 1/T, r, n; u) . 



(4.31) 



Equation (4.29) is a main result of this subsection and 
will be used for the further manipulations. 

It is worthwhile to notice that the functional (4.29) 
has an interesting symmetry. Let us make a shift of the 
variables 



V,^V^-(l-a)^Fr 



(4.32) 



in the fmictional integral (4.29), where a is an arbitrary 
constant. Using Eqs. (4.28), we find for the transforma- 
tion 



2i/V2i J 'ip{X)F {X) dX 2i/V2i J ■ip{X)F (X) dX 
+ {1- a)u j dxdn [h'^{n) + h'^{-n)] , 

(4.33) 

where the notation (2.2) and the convention (4.13) for 
the angular integration are used. Analogously using 

hh = 

that can easily be checked using the definition of h, 
Eq. (3.20), one obtains from Eq. (4.22) 



j V {X) Chip {X) dX^ Jip (X) Chtp {X) dX 



'ijj{X)jC.o'tpCF^ {X)dX. 



(4.34) 



The extra term appearing in Eq. (4.33) has the same 
functional form as the interaction (3.25b) and can be 
incorporated into renormalization of the interaction con- 
stant, whereas the extra term in Eq. (4.34) can be ac- 
commodated into redefinition of the operator F'^{X) 
F^{X;a) in Eq. (4.28b) as 



^|(X;a) = -L(J 



1 



g \ / H \ / eh 

u [adr + i{l — a)vFnV] h^{x, n) 

yU [adr — *(1 — a)vF'n.V] h^{x, ~'^)y 

(4.35) 

Accordingly, the low-energy representation of 
Eq. (3.26c) can be written for an arbitrary param- 
eter a as 

exp (^- ^) = j DhWt ({h} ; a) ({h} ; a) (4.36) 
where [cf. Eqs. (3.25b), (3.25c)] 

Wt(a) = exp | — J h{x,n) f^i(Q;)h (a;, n) rfa;(in| , 



2rt{a) = -f 



1 -I- Q!F<^ ■ 



(4.37) 
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and the angular integration over the whole rf-dimensional 
sphere is meant, and convention (1.2) is implied. 

The partition function Zs{a) is a generalization of 
Eq. (4.29): 



Z,{h,a} = I Dtpexp[-Sf] 

X exp 2iyV2i j^{X)F {X; a) dX 
F{X;a) = F\X) +F^{X;a), 



(4.38) 



where vectors F^;F^(a) are given by Eqs. (4.28a) and 
(4.35), respectively. A particular choice of the parameter 
a is merely a matter of a convenience. 

Closing this subsection, we recast the supersymmetry 
breaking terms in Eq. (4.38) into a form more convenient 
for further application, We introduce a 16-component su- 
pervector J^q 



3^0 = 



V2 



( 



eh 



(4.39) 



and the operator is given by 



i 


(2a- 


1)4 




T++T-; 






(2a- 




- c)' 


T- +T+, 


(4.40) 


i± 












2 













and Co, jCq are defined in Eq. (4.25). The conjugation 

for the matrix operator is given by Eq. (4.20), and the 
supermatrix fa is given by Eq. (4.27). Then it is easy 
to check by exphcit calculation that Eqs. (4.38), (4.28a) 
and (4.35) can be re-written as 

J^{X)F {X; a)dX = J (X) L,„m^{X)J^o dX 

= J To%{X%,afl^^ {X)dX 

(4.41) 



where operator H,. is given by Eq. (4.26), and summation 
over the repeated index "f = x,y,z is implied ™. The 
latter formula is the most convenient for the integration 
over h which will be performed in the next subsection. 



™' When wc write the index in the spin space explicitly, we imply 
that matrix Is should be dropped from definitions (4.27) and 
the relevant supermatrices have the dimensionality 16 x 16 



C. Averaging over the Hubbard-Stratonovich field. 

The argument of the exponential in Eq. (4.38) is a lin- 
ear functional of the Hubbard-Stratonovich field h, and 
thus the integral over h in Eq. (4.36) is purely Gaus- 
sian. The field h enters both the function F and the 
Lagrangian (4.24). This means that the new effective 
field theory will contain quadratic, cubic and quartic in 
xjj terms. The quartic term originates from the averaging 
of the supersymmetric part of the action, and therefore, 
it preserves the supersymmetry, whereas the quadratic 
and cubic terms lift it. 

Performing Gaussian integration over h in Eq. (4.36) 
with the help of Eqs. (4.41) and (4.26), we find the contri- 
bution of the spin modes to the thermodynamic potential 



-Tin 



/ 



exp(-5[V'])£'V 



(4.42) 



with 



S [^] = So {^} + 52 im ; a] + S3 [m ,a]+Si [{^} ; a] . 

(4.43) 

In Eq. (4.43), the free supersymmetric part of the action 
can be written as 



So m = -2iv \ (X) 



Co - iSK 



{X) dX (4.44) 



where Co is given by Eq. (4.25), the summation is implied 
over the repeated spin subscripts 7, [see also footnote af- 
ter Eq. (4.41)], variables X are defined in Eq. (4.23), and 
the convention (4.13) is used for the angular integration. 

The term ^4 describes the quartic interaction and it 
takes the form 

S4 [{ip} ; a] = -Ausspjes/Si-fi ^ Aj^ / dX 

X (V^^ {X) f3n,V7 {X) u) ti (u^p^ {X) fsfl.VTi (^)) 

(4.45) 



where 



A = 



1 1 
1 -1 



(4.46) 



and we introduced the 16 x 16 self-conjugated supermatri- 
ces [sec Eqs. (4.18), (4.27), and footnote after Eq. (4.41)] 

Hi = 1, fl2 = E3, As = Aifs, fi-i = AifsSs. (4.47) 

The significance of matrices II3 4 will become clear in the 
next subsection. 

The operators Fj here are slight modification of Ft (a) 
in Eqs. (4.37), (3.25c): 



fi(Q;) = /7i(a); 



(4.48a) 
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where action of the cut-ofF operator / is defined in 
Eq. (3.25e) and the operators "ii{a) are defined by 

1 

Y)ib]{X)= j dni j dui7j (nrri;M,ui)6(r,T, ni.ui). 



(4.48b) 

Hereinafter, the convention (4.13) is used for the angular 
integration. The kernels in Eq. (4.48b) are given by 



7i (nni;u,Mi;a) 



72 (nni;u,ui;a) = -- 



ni 



a + F" 



(4.49) 



and they are independent on the parameters u,u\. We 
will see that this will change when we consider the fluctu- 
ation corrections to the bare action (4.45). As wc will see 
later, the most interesting effects will come from n ~ rii 
so that the notation 'yf(h) for the forward (backward) 
scattering will be self-evident. 

The tensor Sap-y is the antisymmetric tensor of the 
third rank (£123 = 1) and the summation over repeated 
indices is implied in Eq. (4.45). Also, the relation 



(4.50) 



holds. 

The term [{^} , a] describes the cubic interaction 
and we write it as 

^3 [{-^1 , q\ = -'ii'V2ies0j ^ Xij / dX 

(4.51) 

where operator l^^a and supervector J^q are given by 
Eqs. (4.40) and (4.39) respectively. 

At last, the quadratic term <S2 [{ip} ', a] reads 



S2 m 



2w ^ Xij j 



dX 



X [:Folu,aT3lLji^5 (^)) f , (^Tolu^aTalij^S (^)) • 

(4.52) 

Equations (4.42)- (4.52)) completely specify the field 
theory that describes the collective spin excitations. We 
will see that the interaction between the modes given by 
the terms 53.4, leads in the limit T ^ to logarithmically 
divergent terms of the perturbation theory in these inter- 
actions. These divergencies make the theory non-trivial 
and interesting. The logarithmic contributions can be 
summed up using a renormalization group (RG) theory. 
This will be done in Sec. V. Before doing so, however, we 
will slightly generalize the action to a form reproducing 
itself under the renormalization group procedure. 



D. Further generalization of the theory. 

We start with generalization of the quartic interaction 
by including all the matrices 11^. sec Eq. (4.47), to the 
interaction part of the action. We increase the dimen- 
sionality of matrix A of Eq. (4.46) as 



/I 1 -1 -IX 

1111 
1-1 1 -1 
Vi -1 -1 1 J 



(4.53) 



with the properties 



(4.54) 



Xik = 4(5fci; ^ Xik = 4(5i2- 

i=l fe=l 

Then, Eq. (4.45) can be re-written as 
S4[{ip};a\ = -2uesp-fSsi3i-fi ^ Xij / dX 

X (v^^ (X) f3n,V7 (^) f z (ui^is, (x) f3n, V71 (x)] 

(4.55) 

where the operators Fj are defined by their action on any 
function b{X) as [cf. Eqs. (4.48)] 



Tib 



i 

(X) = Jdm jdui JdrJ{r,) 



X (nni;u,ui;r5^) 6(r -|- ri, r, ni, ui), 

(4.56) 



where 



r_L(n) = r — n (r • n) 



(4.57) 



denotes the coordinate transverse to the direction of 
the momentum, the cut-off function /(r) is defined in 
Eq. (2.29), and [cf. Eq. (4.49)] 



Ti {6; u, ui; r^) = r2 {6; u, m; r^) = j'}i9); 



(4.58) 



Ts {9; u, ui; r^) = T4 {6; u, m; r^) = ^^{9). 
Analogously, we rewrite Eq. (4.51) as 

[{^} ,a] = -2vV2iep^s ^ Xij dX 

i,j = l <T=± •' 

/_ . X , /_ . N (4.59) 

X {X) fgn,^^ {X) uj [D^T^U.^Ps {X)] 

where, similarly to Eq. (4.56) 
Bfb 



(X) = Jdni jdui JdrJin) 


X (nru;u,ui;r^) b{r + ri,T,ni,Ui), 

(4.60) 
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operator lu,a and matrices t± are defined in Eq. (4.40), 
and the bare value of the couphng functions are 

{9; u, ui; r^) = Ti {9; u, m; r^) . (4.61) 

Finally, the quadratic term, Eq. (4.52), is recast as 

X (v^(X)n,r3roa,) Af^^^ (S,,r3fl,V5(^)) , 

(4.62) 

where are defined in Eq. (4.59). (The different overall 
sign in comparison with Eq. (4.52) appears because the 
matrix 73 is anticonjugatc, T3 = —73). 

Though it appears that four coupling matrices 
Af^'^^ i = 1,2,3,4 may be present in Eq. (4.62), only 
two of them actually give a non vanishing contribution. 
Indeed, using Eqs.' (4.59), (4.39) and (4.47), we find 

fisB (g) Sn3 = niD (g) Wtii, fizD (g) Wti2 = 1I4B (g) en4, 

which leads using Eq. (4.53), to 



where P^'^ are supermatrices such as 



i= 1,4. 



^A^feftfeD^Dnfe = 

fc=i 

The non- vanishing operators A^^''^ 
Eq. (4.62) are defined as [cf. Eqs. (4.56), (4.60)] 



(4.63) 

2,3, in 

(4.64) 



X A'^'"^ {nni;u,ui;rj^) 6(r + ri, r, ni, Ui). 

Equations (4.62)-(4.64) reproduce Eq. (4.52) for 

A^i"= {9; u, ui; r^) = Ti {9; u, r^) , 1 = 2, 3; 

(4.65) 

We will see that the form of the action (4.55), (4.59), 
(4.62) for q; = 1/2 will be reproduced by the renormal- 
ization group but the relation between constants (4.58), 
(4.61), (4.65) will be violated, so that Eqs. (4.58), (4.61), 
(4.65) will serve as initial conditions for the renormaliza- 
tion group fiow of Sec. V. 

The reason for the introduction of the 11^ matrices 
(4.47) into the definition of the interaction actions is that 
they separate the combination of the supervectors which 
may transform the free Lagrangian £0, Eq. (4.25), into 
t\. We will in the next section that it will be a necessary 
condition to give rise to the logarithmic divergence. 

To understand such partition, note that any superma- 
trix P can be represented as 



(4.66) 



= 0, {p(i\f3Ai}=0 (4.67) 

= 0, [P^^\T:ih] =0 

{P(^),S3}=0, {p(^),r3A,}=0 

where [...,...] stands for the commutator, {...,...} for 

the anticommutator, and the relevant supermatrices are 
defined in Eq. (4.27). 

It is not difficult to invert Eq. (4.66): 



1 4 



(4.68) 



k=l 



where the 4x4 matrix A is given by Eq. (4.53). Eq. 
(4.68) can easily be checked by using the property (4.54) 
and the commutation relations of the matrices VLk- 
One can see that 

Str (pWpO)) = SijStv (pW)' , 

Str (pq) = Str (pWg^'^)) (4.69) 



for arbitrary supermatrices P and (j, where the super- 
trace operation is defined^'^ as 



Str 



a p 
a b 



= Tra- Tr6. 



(4.70) 



i=l 



We note that the following very useful relation 

4 

Y aii&i2^nfei^i2fe2Str (^Atlk^tlk^Btlk^ilk^'j 

ii,i2,fei,fe2 = l 

4 

= 4 ^ ai6,A,feStr (^AUkBUk) (4.71) 

i,fe=l 

is valid for arbitrary coefficients Ci and bi. One can prove 

Eq. (4.71) by a direct calculation using the fact that a 
product of two matrices (4.47) is once again one of the 
matrices (4.47). Finally, combining relations (4.68) and 
(4.71), one finds 

e 4 

str (p(^)Q(^)) = ^J2 ^i'^^*^ (P^kQ^k) (4.72) 

k=l 



V. PERTURBATION THEORY AND 
RENORMALIZATION GROUP. 

This section contains the perturbative analysis of the 
field theory derived in the previous section. We will start 
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in Sec. VA with the brief formulation of the rules of 
the diagrammatic technique^^ emphasizing the aspects 
different from the conventional models. We will show 
the origin of the logarithmic divergence in Sec. VB and 
formulate the renormalization group (RG) procedure of 
summation the leading logarithmic series in Sec. V C. We 
will be able to derive the RG equation for the coupling 
constant Eqs. (4.55), (4.59), (4.62). The solution of the 
RG equation will be done in Sec. VI. 



A. Rules of perturbation theory. 

As usual, we would like to construct the expansion 
of the observable quantities in terms of the interaction 
vertices produced by Eqs. (4.55), (4.59), (4.62) and the 
Green functions (in our case supermatrices) of the free 
motion 



Let us discuss peculiar features of the perturbation the- 
ory. First of all, due to the supersymmetry of Sq, the 
averages of the operators Ai not perturbing the super- 
symmetry vanish 



0; 



0; 



(5.4) 

This leads to the cancellation of the closed loop contri- 
butions, see e.g. Fig. 6b). 

Second feature originates from the dependence of the 
supervectors -0, -0 on each other, see Eq. (4.19). This 
mutual dependence makes the rules of the Wick contrac- 
tions of the supervectors very similar to ones for the real 
fields (so that the arrow in the Green function looses its 
meaning). To illustrate this point, consider the connected 
average (Qi,2 are the arbitrary self- conjugated superma- 
trices, (5l,2 = <9l,2 )'■ 



ao(Xi,X2) = -4i!/(V(Xi)«)V(X2))o (5.1) (V02^))^ = (^Oi'«/')(V'Q2V') (5.5a) 



where variables X are defined in Eq. (4.23), and the av- 
eraging means 



(•••)o 



. . . exp {-So 



(5.2) 



the free action is defined in Eq. (4.44), and the normal- 
ization is trivial due to the supersymmetry: 



+ {tPQiiP){tPQ2il^), (5.5b) 



where over- and underbrackets stand for the Wick con- 
tractions. The line (5.5b) can be transformed with the 
help of Eq. (4.21) as 



/ 



exp(-5o[V'])£'V' = 1- 



The factor — 4ii^ is introduced in Eq. (5.2) for the sake of 
convenience. All the higher order averages are, then, to 
be found using the Wick theorem and Eq. (5.1). 
Using Eqs. (4.44), (4.25), one easily finds 



Go {Xi,X2) =S{ui- U2)S (ni;n2) 

x^E / ^e--"(--^^)+'''('-— -=)Go(a;„,k;ni); 

Go (w, k; n) = 



icoAi + UirknfsSa -|- i6A 

a;2 + t;2(k„)2+^2 • 



(5.3) 



From Eq. (5.3) we see that terms involving 6^0 are dan- 
gerous only for zero Matsubara frequency contribution. 
These contributions are associated with the real scatter- 
ing event with the energy transfer much smaller than 
temperature. Such processes, though determining the 
kinetic of the system, are not interesting for equilibrium 
thermodynamics, and will be considered elsewhere^''. For 
the logarithmic contributions considered further in this 
paper, the real processes are not important and we will 
put 6 = from now on. 



(5.56)= (tAOit/>)O0Q2t/'), 

which coincides with the first term in the right-hand-side 
of Eq. (5.5), because Qi is self-conjugate. As the result, 
we obtain 

^ (V^Qi V) (?(32V) )^ = -2 Str [QiaoQ2ao] , (5.6) 

where we omitted the trivial factors of proportionality 
between the averages and the Green fimction, Eq. (5.1). 
The appearance of the factor of 2 in such an average is 
the feature of the real fields. The minus sign in Eq. (5.6) 

originates from the definition of the supertrace operation, 
Eq. (4.70), where the commuting sector is taken with the 
negative sign. 

To further utilize the analogy with the real fields, let 
us consider a connected average involving eight fields but 
with only four fields contracted (such an averaging ap- 
pears e.g. as a correction to the interaction constant for 
the quartic scattering term) 



Ji = 

J2 = 



(t/.gitA)(^Oit/>) {^PQ2 1P){^l>Q2'^P) 



(5.7) 



If were usual fermionic field, contributions J^i and J^2 
would be responsible for different processes (particle-hole 
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and particle-particle, respectively). For the superfields, 
however, we can transform Ji [cf. derivation of Eq. (5.6)] 
and thus obtain 



;. l< . II, u 



i.e. the construction describes precisely the same con- 
tributions. Analogously one finds 



(l/>QlV)(V'On/') (^/'Q2V')(V'02V') 



so that the permutations of the fields in the vertices do 
not give rise to any new effects but simply lead to the 
multiplication by a factor of 4 - number of trivial sym- 
metries of the interaction vertex. That decrease of the 
number of different contractions is a great simplification 
in a further derivation. 

This observation enable us to formulate the simple di- 
agrammatic rules for the generation of the perturbation 
expansion, see Fig. 4. 

As usual the summation of the non-fixed by the ex- 
ternal legs or conservation laws coordinates oj^, k^, n^, u^, 
must be performed, and the supermatrix product to be 
calculated. The closed loop of the supersymmetric Green 
functions bring — Str of the product of all the terms in 
the loop, [cf. comment after Eq. (5.6)]. 



B. Identification of logarithmic divergence. 

Having established the basic rules of the diagrammatic 
technique we are ready to demonstrate the logarithmic 
singularity in any dimensions. 

As usual in the logarithmic series, one has to look at 
the perturbation theory for the interaction vertices. The 
lowest order diagram of interest is shown on Fig. 5 and 
all the remaining terms are enumerated and discussed in 
Sec. VC. 

According to the rules of the diagrammatic technique. 
Fig. 4, the term of Fig. 5, (let us denote it by jTs) can 
schematically be presented as 



AT 

V 



k 

X 7i(nrn5)7j(fun5) 

X Str(f3(5iGo (<^n, k; ni) TSQ2G0 {ojn, k; n2)), 

(5.8) 

where the supermatrices Qi,2 indicate all the combina- 
tions of the supermatrices and direct products of the su- 
pervectors standing on the left and on the right of the 
Green fimctions, see Fig. 5. The only important point is 
that we will neglect their momentum dependence, which 



£Go (w,k, n) 



7 

Ml 



71, k 

/3i,k-hq 

4 



Ml 

k 



5 

q 



i,j=l a=± / 



M,k,n,7 




Mi,k, ni,7 



CT1,2=± 



FIG. 4: Basic element for the diagrammatic technique, Ex- 
ternal legs of the vertex arc amputated. The expression for 
the interaction vertices are found from Eqs. (4.55), (4.59), 
(4.62), and the additional symmetry factors (in brackets) 
are discussed in text. The lines in the analytic expressions 
for the diagrams indicate the directions of the insertion of 
the Green function Go- Expressions for D± are obtained 
from Eqs. (4.59) and (4.40) with £0 = iwAi -|- VF(kn)E3-r3, 



-iwAi -h VF(kn)E3f3. Cut-off function / is defined in 



Eq. (2.13). 



suffices our aim for the logarithmic accuracy. The ex- 
tra matrix fs and the numerical coefficient in front are 
introduced for the convenience. 

Substituting Eq. (5.3) into Eq. (5.8), and keeping only 
non-vanishing terms, we find 



7,(nln5)7j((nln^)) 



/(k-k/)/(k-k, 



X (4fc|fc|strQiS3Q2S3 - w2StrQif3Aif3g2Ai) ; 
2 = k • ni,2. 



(5.9) 



23 



kf, ojf,ui,ni \ k,a),Mi,ni / k,, a;i,Mi,ni 



k/, w/,W2,n2 




FIG. 5: Lowest order correction to the quartic interaction 

vertex (external legs arc amputated). This contribution log- 
arithmically diverge at T ^ 0,ni ^ n2, independently on 
dimensionality. Objects separated by dotted arcs are named 
Qi,2 in Eq. (5.8). 



Direct examination of Eq. (5.9) shows that the integral 
can exhibit logarithmic divergence for ni n2 only if the 
two terms in the last factor have the same sign. The last 
condition is conveniently taken care of by representing 
each Q matrix in the form of Eq.(4.66), and using the 
relations (4.67). To facilitate further manipulations, we 
introduce 



n=(ni+n2)/2; (5n = ni - n2; 
fell = k • n; k_L = k — fciin, 



and consider \6n\ <C 1. After integration over in 
Eq. (5.9), one finds [using the decomposition (4.66)] 



J5*^=7i(|H)7i(l<5n|)Str Qp^Qf 



d''-iki4u;„/(ki-k^)/(k^-k,) 



(2vr)" 



(5.11) 



where the upper limit on the frequency summation ap- 
pears because we neglected fc|| dependence of the cut-off 
function / and the typical momenta contributing to the 
integral are of the order of uin/vp- 

Summation over Matsubara frequency leads immedi- 
ately to a logarithmic result independently on the di- 
mensionality of the system: 



Qi ''Q2 



x/(2)(k/;k01n 



roe 



(5.12) 



where the infrared cut-off of the logarithm is determined 

by 



i = max {T, \Sn\ vfTq ^} <C 



Vf_ 

ro ' 



(5.13) 



The supertrace in Eq. (5.12) can be re-expressed using 
Eq. (4.72) as 



Str 



(3) 
2 



Qi 



-lEE(-l)'^^^'^Str 



i=l,3 fc=l 



QiUkQzUk 



(5.14) 



Although the logarithmic divergence (5.12) is present 
in any dimension, the coefficient in front of the logarithm 
in Eq. (5.12) is determined by the dimensionality and the 
details of the ultraviolet cutoff: 



^/(k^ - kj)/(k^ - k,). 



(5.15) 

where the cut-off function /(k) is given by Eq. (2.13). 
For d = 1 there is no integration over the transverse 
momentum, f{k) cuts the logarithmic divergence only, 
so 



f(2) 
Jd=l 



Ml 



(5.16a) 



for k/^i < (l/''o) and decreases rapidly for the larger 
momenta. For c? = 2,3 we notice that logarithmic con- 
tributions originate from the region |fc|[| <C |k_L| <C 1/ro, 
and this feature will persist in all the further terms of the 
perturbation theory. Neglecting the parallel components 
in Eq. (5.15), we find 



(5.10) /^'^(k/;k,)=Mdy 



,i(ki-kj)r 



112 = 4: {pFroy^ ; Ms = 47r (pFro)"^ 
ro 



fi 



f- 



(27r)'^-i 



/(k), 



(5.16b) 



where the coordinate integration is in the plane r • n = 0, 
The significance of f± is the regularization of the fields 
that are otherwise singular functions of the transverse 
coordinates. 

Equation (5.12) demonstrates that the field theory un- 
der study is logarithmic in any dimensions. The correc- 
tions coming from the interaction diverge in the limit 
{T, \ dn\} 0. Therefore, we can use a renormalization 
group scheme for calculation of physical quantities. This 
can be done in the limit of small F considered here. 

A remarkable feature of the corrections is that the main 
contribution comes from configurations with either par- 
allel or antiparallel alignment of the vectors n. To some 
extent, the spin degrees of freedom of the electron system 
have a tendency to forming a one dimensional structure 
and this happens in all dimensions. 



C. Integration over fast variables 

This subsection is devoted to summation of the per- 
turbation series in the leading logarithmic approxima- 
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tion. It means that the expansion for a physical quan- 
tity y is classified not in powers of a coupling con- 
stant y — Yin ^"7" but as an expansion of the type 
y = X]„ [7ln(. . . )]" a„(7). The renormalization group 
(RG) corresponds to the Taylor series expansion of each 
function a„ (7) , whereas the value of the logarithmic fac- 
tor itself can be large. 

Following the conventional scheme, see e.g., Ref. 58. 
wc subdivide the supervectors ij} into a slow, ^ (X), and 
fast, t {X), parts 



il,{X) = ^ {X) + T (X) , 



(5.17) 



and integrate ovca- the fast variable T {X) using the per- 
turbation theory in the effective interaction. 

One should be careful defining the fast and slow vari- 
ables because, as we saw from Eqs. (5.9) (5.12), the loga- 
rithmic contributions originate from the configuration of 
the fields highly anisotropic in space (smooth along di- 
rections of the momentum n and sharp in the transverse 
direction). Therefore, we can not separate the fast and 
slow variables considering the moduli of the momenta'^. 

Fortunately, this problem can be avoided because we 
can define the fast and slow variables with respect to 
frequencies only. As the main contribution in the integral 
over fc||, see Eqs. (5.9), (5.11), comes from k\\ ~ i^j/vp, 
this type of the separation is suSicient. As concerns the 
perpendicular components fcj^ , they do not participate in 
the renormalization group treatment entering equations 
as a parameter (like the other variable u). 

After decomposition (5.17) the action acquires the 
form 

S {*, T} = <So {*} + {T} + S2 {V-} + <Sint {*, T} , 

(5.18) 

where 

5i„t T} = S2 {T} + 53 {* + T} + Si {* + Y} . 

(5.19) 

The free action for the fast fields has the form, cf. 
Eq. (4.44), 



{T} = -2iv 



T-y C-o — i {xoJc) ^ T^dX, 



(5.20) 



where Co is given by Eq. (4.25). The second term in 
brackets leaves the contribution only from frequencies 



< 



(5.21) 



where uic is the running cutoff'-of the problem, so that 
Iwl < iVc and x < 1. 



With such a choice, one step of renormalization group 
transforms the running cut-off as 



(5.22) 



Our goal is to obtain the correction to the action of 
the slow variables \E' arising due to the interaction with 
the fast fluctuations: 

6S^ = - In (exp [-5int T}])^ - 5int 0} . (5.23) 

Hereinafter, the averaging over the fast fields T is defined 
as cf. Eq. (5.2) 



{■■■h 



J ... exp{-So{r})Dr. 



(5.24) 



The integration over the fast field {T} is performed 

using the Wick theorem and, thus, all the machinery of 
VA is still applicable. The only difference is that in the 
intermediate lines one has to replace Go Gr, where 



Gr (w, k; n) 



(5.25) 



which differs from Eq. (5.25) by the regularization term 
restricting the domain of the frequency integration from 

below, Awc, and from above /i (^i^)- Smooth function 

fi{x) has the asymptotic behavior /i(a;) — > 1, a; <C 1 and 
f{x^ 00) ^0"". 

Analogously to the bare action, its correction can be 
decomposed 

SS^ = SS4 + 6S3 + 6S2 + 5Sq. (5.26) 
We will consider each of those contributions separately. 

1. Renormalization of the quartic term, dS4 {^} ■ 

The first loop diagrams leading to the renormalization 
of the quartic interaction are shown in Fig. 6. Only the 
diagram. Fig. 6a), may produce a logarithm. Indeed, the 
diagram Fig. 6b) contains a closed loop and vanishes be- 
cause of supersymmetry, see Eq. (5.4). The diagram 6c) 
does not produce a logarithm because of the locations of 
the poles in the corresponding Green functions, as it will 
be more formally discussed in the end of this subsubsec- 
tion. 

To obtain the analytic expression for the diagram 6a), 
we apply the rules of Fig. 4 and notice that the result 



It should be contrasted to the consideration of Ref. 5 where RG 

language was merely used to reformulate known^'-'^'''^^ results. 
In the RG language, the Fermi liquid constants arise in d > 1 as 
the dimcnsionally irrelevant couplings and, thus, RG procedure 
for calculating their values is useless. 



^ The particular functional form of the cut-off will be not im- 
portant in at least first loop renormalization group calculation, 
though it might be that the more accurate choice will be required 
for the higher loop calculations [we did not investigate such loops 
in details]. 
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oj, k, n, u 



-= ^GT(a;,k,n) 



= -<S4 
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b) 



m 



FIG. 6: First loop renormalization of the quartic interaction. 

The notation is introduced in Fig. 4, and filling moans the ver- 
tices with rcnorrnalizcd interaction constants 7^. The Green 
function for the fast fields, Gr, is introduced in Eq. (5.25). 
Factor [1 / 4] in diagram a) account for the symmetries of the 
diagram. Diagram (b) vanishes because of the supersymme- 
try, see Eq. (5.4). Diagram c) does not produce logarithmic 
contribution, see Eqs. (5.31) - (5.30). 



has the structure of Eq. (5.8), see also Fig. 5, with the 
replacement of Go Gr and 



1 * , ^ 

Ql,2 = \fv / ^ ^ikeal3i5i^a02S2 



(5.27) 



we find 

4 » 

SSi = 2v^fX ^ \, / dX (*^, {X)%ti,^p, {X)u) 



(5.28) 

where the action of the operator 6Ti is defined by 
Eq. (4.56) with the kernels 



= STi = 0; 
STi (6';M,Mi;r_L) 



-f^duuif^ (r±) In K [Ti (. 



{0; u, ui; r±) = iJ.dUUif±{r±) In x [Ei (. . . )]^ , 

(5.29) 

where we suppressed the arguments in the right hand side 
implying that they arc the same as in the left hand side. 
Equation (5.29) is valid for \9\ < iVcro/vp, otherwise the 
logarithmic renormalization vanish. The function /_L(rj^) 
is defined in Eq. (5.16b). 

The tensor 5^^^^ is given by 



_ 9, 

"/3i7i ~ ^tai/3ii5itai/92<52ta2/3i<5ita2/32i52 
= 2(5-^172^/31/32 + 2(^71/31(572/32 
= ea/3i7iea/32 72 + ""/3i7i ' 
^^'kVi ~ 2(57i/3i (572/32 + ^7i72^/3i/32 + '57i/32'5/3i72 • 



(5.30) 



As matrices (4.47) are self-conjugate and 73 is anticonju- 
gate, one finds using Eqs. (4.20)- (4.21) 

(*/3i (X)f3n,M/7i (X)) = - (*7i (X)f3n,*/3i (X)) , 

and the contribution proportional to S'E will vanish after 
substitution in Eq. (5.31). 

Therefore, the resulting action (5.28) is nothing but 
the original quartic interaction (4.55) with the couplings 
Ti renormalized according to Eq. (4.65). This is sufficient 
to write down the renormalization group equation. We 
will do it in the next section after we consider the trans- 
formation of the remaining terms in the action under the 
RG step. 

Closing our consideration of the quartic interaction, let 
us give a formal proof that the diagram Fig. 6c does not 
give a logarithmic contribution. Once again, we apply 
the rules of Fig. 4 and notice that the result has the 
structure of Eq. (5.8) with [cf. Eq. (5.27)] 



where the indices in the spin space are written explicitly, 
and index i = 1,2,3,4 labels the coupling constant in 
Eq. (5.8). 

It is easy to see that the appearance of the cut-offs in 

Eq. (5.25) leads to the replacement of In — > ln(l/x) 

in Eq. (5.12) without affecting the matrix structure of the 
latter. Using Eq. (5.14) and applying Eq. (4.71) twice. 



Q2 = \fvu\ Aifeea/3i5iea/32 52 

fc=l 

X nfc*/3i ® ¥/32nfef3; 

4 

Qi _ = \fvu\ Aifcec,/3ijiea/32fe 
fc=i 

X fife (*^2nfe'r3*/3i 



(5.31) 
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Using Eq. (5.14), we find 

4 

i=l,3fe=l 

4 

i=l,3fe=l 



nig2] = 0, 



(5.32) 



as the matrices Ilk, see Eq. (4.47), commute with each 

1 2 



other, 
used. 



1, and second of the properties (4.54) is 



2. Renormalization of the cubic term, SS3 {^'} . 

The first loop diagrams leading to the renormalization 
of the cubic interaction are shown on Fig. 7. Similarly 
to what we saw when calculating SS^. only the diagram 
Fig. 7a) produces a logarithm and we turn to the calcu- 
lation of this contribution now. 



= -S3 






a; + n, q + k 




c) 




d) 



FIG. 7: First loop renormalization of the cubic interaction. 
The notation is introduced in Figs. 4, 6 and the filling means 
the vertices with renormalized interaction constants Pi. Di- 
agram (b) vanishes due to the supcrsymmotry, similarly to 
Fig. 6b. Diagrams c) d) do not contain logarithmic contribu- 
tion, similarly to Fig. 6c. 



the quartic term is the presence of the differential opera- 
tor ID- in the expression for the vertex. Using Eqs. (4.59) 
and (4.40) we write 



u- 



(2a - 1) £o(w, k) - Cl{oj -h 1^, q -h k) 



--^o£S(-f2,-k)f_ 



(5.33) 



where frequencies and momenta are arranged as in 
Fig. 7a). 

The terms in the second line cancel the small denomi- 
nator in one of the Green function. Integration over uj, k 
does not produce terms cx In >ir because 



<kvdk]\ABG-r{ijj,k\\) ~ — i (kwc In k) A, 



(5.34) 



i.e, it is determined by the lower limit of the integration 
and must be excluded from RG scheme. 

The term in the last line is not affected by the; integra- 
tion, so the result can be once again recast in the form 
of Eq. (5.8), and the calculation proceeds similarly as it 
was done for the quartic term. Instead of Eq. (5.27), we 
find 



fe=l 



k=l 

where we introduced the notation 

D,, =!<,(« = 1/2) . 



(5.35) 



(5.36) 



Repeating the same steps as when deriving Eq. (5.28) 
and using the identity 

(ea/3i5i£^a52/32) ^i52 725i = ^/Si&T) (5.37) 

we obtain, [cf. Eq. (4.59)] 

5St = -2iyV2ief3^s E E / 

X (^^0 (X) fsftj* ^ (X) u) SBt [W+Tsflj^s (X)) , 
6S^ = 4uV2ief3.ys ^ij / dXu (5.38c) 

i,7 = l a=±'' 



(5.38a) 
(5.38b) 



The only difference in this calculation from the one for ^ 



n,*^ {X) ^ dB-B_f3^ \Uj^s {X) (X) fa 
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where the sign difference between Eqs. (5.38b) and 
(5.38c) appears because of the definition of the super- 
trace (4.70). 

The action of the operators 5Bf in Eqs. (5.38) arc 
defined by Eq. (4.60) with the kernels given by [cf. 
Eq. (5.29)] 



SBf = SBf = 0; 

SB+ = -2^^dUUlh{r^) Inx TiB+; 
SB^ = -lJ.dUUif±{r±) Inx TiBf 
SB^ = 2/idUUi/_L(r_L) In x TsB^; 
SB^ = HdUUif±{r±)\nM: T^B^ 



(5.39) 



for \6\ < LOcro/vp, otherwise the logarithmic rcnormal- 
izations vanish. We did not write the arguments of the 
kernels implying that they are the same as in the left- 
hand side of Eq. (5.29). 

Equation (5.38c) is apparently not of the original form 
(4.59) yet as the differentiation in Eq. (5.38c) acts on two 
fields on its right whereas the derivative of only one field 
is present in Eq. (4.59). However, it can be transformed 
using 



(as matrix fs is anticonjugated), and z denote either r 
or r. We, thus, re-write (5.38c) as 



5So = 4iyV2ii 



dX 



(^0 {X) faft,*^ (X) u) SBr (b.tsIIj^s {X)) 



InUj^s {X) 

(5.40) 



where the notation ED means that the differential opera- 
tors included in P act on the left: 



cd zd= [dzc) d 



(5.41) 



for arbitrary functions c, d. After integration by parts in 
the last term we obtain 



5So = 2vV2ii 



4 



dX 



X (X) {X) 5BI (n.T^Iij^s {X) 

(5.42) 

that has the same form as Eq. (4.59) for the particular 
choice of the parameter a 



1 



a = 



(5.43) 



b) 



/ 



\ 



'///////, 



FIG. 8: First order corrections to the quadratic term. Di- 
agram a) vanishes due to the supcrsymmotry. Diagram b) 
vanishes after w,k integration, see Eq. (5.34). 



In what follows wc will use only this value of a. 

Equations (5.38), (5.42) shows that the cubic term is 
reproduced under the RG step and Eq. (5.39) determines 
the new values of the coupling constants. 



3. Renormalization of the quadratic interaction, 6S2 {^} 
and free action SSo {^'} . 

The one loop diagrams that may change the values of 
^50.2 {^} arc shown on Figs. 8 10. 

One immediately notices that all the diagrams of the 
first order in quartic interactions. Fig. 8, vanish. Indeed, 
diagram of Fig. 8a) vanishes because of the supcrsym- 
metry. Diagram of Fig. 8b) vanishes because it involves 
integration of one Green function only and cannot pro- 
duce a logarithmic divergence, see Eq. (5.34). 

Among the diagrams of the second order. Figs. 9,10, 
only Fig. 9 a,b) gives the logarithmic contribution into 
the quadratic interaction. We transform contributions 
involving D_(cj,k), I}-{u>,'k) on diagrams Fig. 9a ac- 
cording to Eq. (5.33), for a = 1/2. Then, the contribu- 
tions from the last line of Eq. (5.33) are not affected by 
the integration, so the result can be once again recast in 
the form of Eq. (5.8) with [cf. Eqs. (5.35), (5.27)] 



Qi 



^ fc=i 

(5.44) 

The first line of Eq. (5.33) for the quadratic in- 
teraction does not produce any logarithmic diver- 
gence for terms involving either I!'+(cl', k) . . . ID)_ (w, k) or 
B+(w,k) . . .D_(_a;,k) due to the integral Eq. (5.34). The 
term involving D_ (w + fi, k + q) . . . D_ [yj. k) have all the 
denominators cancelled and the integral is ultraviolet di- 
vergent. This ultraviolet divergence will be discussed 



Other choices of a. would require additional shifts of the fields 
in order to reproduce the cubic term in the RG. Therefore, the 
choice (5.43) is the most convenient for the sake of the calculation 
though the final physical answers can not depend on a. at all 
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= -<S2 



a) 




c) 
+ 



e) 




FIG. 9: First loop corrections to the quadratic term <S2. The 
notation is introduced in Figs. 4, 6 and the filhng means the 
vertices with renormahzed interaction constants /S.l^"'^ . Dia- 
grams a,b) arc logarithmic. Diagram c) vanishes due to the 
supersymmetry, similarly to Figs. 6b, 7b. Diagrams d) e) do 
not contain the logarithmic contribution, similarly to Figs. 6c, 
7c,d. The factor [1/2] on the panels a,b) appears because of 
the symmetry of the diagrams, 



later on in this subsection after we complete the deriva- 
tion of the logarithmic terms. 



In the diagram Fig. 9b only the terms involving 
D-(-(a;, k) . . . D+(w, k) are logarithmic. Diagrams involv- 
ing 5+ [u, k) . . . ID_ (o), k) and D_ [uj, k) . . . 1D)+ (w, k) are 
of the type of Eq. (5.34) and do not contribute to the 
logarithmic renormalization group. Finally, the term 
D_(a;,k) . . .D_((jj,k) gives rise to the ultraviolet diver- 
gence which will be discussed shortly, see Eq. (5.51). The 



logarithmic part of Fig. 9b is of the form (5.8) with 

4 



Q2 



5i52 



A 

1 



fe=i 



(5.45) 



Collecting all the logarithmic contributions from 
Fig. 9a,b) with the help of Eqs. (5.8), (5.44), (5.45) and 
the identity 



we obtain [cf. Eq. (4.62)] 

4 



(5.46) 



6S2 [{i^} ; a] = -w Xij J 

i,j=l 0'l,2==t 



dX 



(5.47) 

The action of the operators S^f^ for i = 2, 3 are given 
by Eq. (4.64) with the kernels [cf. Eqs. (5.29), (5.39)] 



SA^'"^ = 0; 

5A^+ = 2jUdUUi/j_(r_L)lnx 



r3A++ + (B3+)' 



(5.48) 



SA+- 



and 



2/idUUi/_L(r_L)lnx B3 B^. 



SA^ = 2^dMMi/j.(r_L) ln>< (S3 )^ . 



(5.49) 



Equation (5.48) is valid for \0\ < LOcro/vp, otherwise, the 
logarithmic renormalizations vanish. Once again, we did 
not write the arguments of the kernels implying that they 
are the same as those in the left-hand side of Eq. (5.29). 
The reason why the correction ^A3 ~ is written sepa- 
rately from all the other couplings will be explained mo- 
mentarily. 

To complete the calculation of the correction to the 

quadratic interaction, wc have to compute actually the 
ultraviolet divergent terms in 9a-b. We found 

[Fig. 9a + Fig. 96]„„ ^ ^ Oii^l) { [/3r] ^ - ^.A" } . 

i=2,3 

(5.50) 

The ultraviolet divergences cancel each other for the ini- 
tial couplings (4.61), (4.65). In fact, the vanishing of 
such divergences precludes the formation of the gap in 
the spectrum of the excitations forbidden by the spin ro- 
tational symmetry, and should be valid in any order of 
the perturbation theory. On the other hand, the accuracy 
of our renormalization group procedure does not allow us 
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ni 




n, M n, w n,u 



FIG. 10: First loop correction to the free action, Sq. This 
correction is logarithmic only for d = 1. For d > 1, the 
logarithmic divergence vanish due to the integration over rml • 



reproduced after integration over the fast variables and, 
moreover, describe the changes of the of the coupling con- 
stants in Eqs. (4.55), (4.59), (4.62), under the renormal- 
ization. This will enable us to write proper renormaliza- 
tion group equations in a standard way. These equations 
and their solutions are presented in the next Section. 



to the determine finite logartithmic terms from the un- 
certainty "oo — oo". This makes the correction (5.49) 
meaningless. Fortunately, we can use the symmetry of 
the system forbidding the formation of such ultraviolet 
divergences to fix this coupling constant. Requiring the 
most divergent part of the ultraviolet divergence to can- 
cel at any stage of the RG procedure, we find 



r..A; 



(5.51) 



and use this equation for the further RG flow. 

The last interesting diagram is shown in Fig. 10. As it 

conserves both the momentum direction n and the value 
of the coordinates w, k, m, it is natural to classify it as the 
correction to the free action, SSq. The pole structure of 
the Green functions allows for the logarithmic divergence 
at rii n. Calculation is performed using the formula 



(b^^ (w, k, n; u)IlkQo{'^, k, n)Po.2 K n; uj^ 
K++ = K~ = 0; 



(a;,k, n;u) = "'"(a;,k, n; 



u) = u 



Vpin ■ k) 



(w, k, n; w) = K2 ^(w, k, n; u) = w 



+ Vp{n ■ k)2 ' 
2ivFU!{n ■ k) 
a;2 + v%{n ■ k)2 ' 
(5.52) 



that can be checked directly using definitions (5.3) and 
(4.59) for a = 1/2. However, for d = 2,3 integration over 
ni cancels out this logarithmic divergence. For d= 1, one 
finds 



J *7 (X) 



CqR 



^-y {X) dX, 



R ~ B+i3f 1 - Ha + B+B^ 1 + H; 



This correction can be eliminated by the rescaling of the 
fields 'S' ^ [1 — {In >c)R]'^. This rescaling will give the 
third order correction to the coupling constants in the 
interaction part of the action, and thus has to be taken 
into account only in the two loop RG equation. As we 
do not consider such loop in the present paper, we will 
have to disregard SSq even for d = 1. 

Equations (5.28), (4.65), (5.38) - (5.42), (5.47)- (5.48) 
are the main results of this section. They show that the 
functional form of the interaction part of the action is 



VI. RENORMALIZATION GROUP EQUATIONS 
AND THEIR SOLUTION. 



A. General structure of RG equations. 

We have demonstrated in the previous section that 
the functional form of the interaction part of the action 
(4.55), (4.59), (4.62), is reproduced after integration over 
the fast variables, T. On the other hand, each integra- 
tion over the fast variables corresponds to the transform 
(5.22) of the high-energy cut-off 



In Wc — > In >f + In ujc- 



(6.1) 



This enables us to write a most general form of the RG 
equations as 



dV^ 
din LOc 
_dA_ 
dlnwc 
_dAi_ 
dlnuJr 



= 5Ba, (f,;4;4)- 



(6.2) 



The renormalization group fiow starts from ~ vf/tq 
with the initial conditions (4.58), (4.61), (4.65) and it 
should stop at llJc — max(r, \6\vp/r{)). One sees imme- 
diately a significant difference between the standard RG 
scheme and the problem in hand. In our case, the entire 
coupling operators may be renormalized, its renormaliza- 
tion is a functional of all the other operators, etc. Thus, 
we are dealing with the functional renormalization group. 

Surprisingly, in the one loop approximation, the func- 
tional RG equations can be obtained explicitly and solved 
in a closed form. 



B. One loop RG equations. 



The one loop equations for the kernels Fj, Bi, Aj in 
Eqs. (4.56), (4.60), (4.64) are obtained by dividing both 

sides of Eqs. (5.29), (5.39), (5.48), and (5.51) by In and 
taking the limit In = d In ^ 0. As the result, we find 
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(4.58). Comparing Eq. (6.4) with Eqs. (6.3), we obtain 



a In tOc 
dBt{...) 



d In Lu'c 

dr^ dBl 



dhuoc d In 



= -2/.dW±(r±)ri (...)B+ (...); 

= -/idMWi/j_(r±)ri (...) (...) ; 
= 0. 



(6.3a) 



for the set of couplings not having counterpart in the 
quadratic part of the action. Here (...) is the short hand 
notation for the omitted arguments {9]u,Ui;y±). 

For the coupUngs affecting the quadratic part of the 
action, we find 

— = MdUUi/j-ir-L) (6';M,ui;r_L)J ; 

a In LOc 



dBt{...) 
dlnwc 

dlnwc 



= 2/Xdum/±(r±)r3(...)B| (...); 
= /XdWWi/j_(r_L)r3 (...) (...) ; 



dAt+ (...) 



2/idUUi/j_(r_L) 



dXnuJc 

x{r3(...)A++(...)+[B3+(.. .)]'}; 

rfA3-+(...) _ rfA3+-(...) 

dlnujc d In Wc 

= 2^JiauuJ^{v^)Bt{...)B^ {...); 

Ts (...)A3--(...)=[S3-(...)]^ 



rflnti 



dB^ dA?i''= 



rflnt 



0. 



(6.3b) 



The form of Eqs. (6.3) suggests immediately the fol- 
lowing scaling form for the coupling kernels 

Vi {9; u, ui; r^) = 7i {9; u, ui; r^) ; 7°(^)] 

{0; u, uv, r±) = {9; u, uv,v±) ; 'y^{9)] 
Af"^ {9; u, ui; r^) = A^^-^^ (^; w, ur,ri_) ; 7°(^)] 



^ (6*; u, ui; r_L) = uui/Ud/±(r_L) In 



mm 



1 VF \ 

9'roTj 



(6.4) 



dji (0 



dri4__dM 
d^ ~ d^ 



= [7i(0]'; 

= 271 (0/3i+(0; 

= 71 (0/3r(e); 

= 0, 



(6.5a) 



and 



^73 (0 

^/33+ (0 



= -[73(0]'; 

= -273(0/33+(0; 

= -73(e)/?3"(e); 



M = _2A3..(073(0-2[/^3-(0]^ ^^-^^^ 



dA3-+ (0 _ rfA+- (C) 



2/33"(0/?3+(0 



A3-- (073(0 = [/53-(e)]- 



d£, d£, d£, 



0, 



where 7° = 7° = 7/, 7° = 7° = 7;,, see Eqs. (4.49), 



Equations (6.5a) - (6.5b) have to be solved with the ini- 
tial conditions [cf. Eqs. (4.58), (4.61), (4.65)] 

7i(^ = 0) = (e = 0) = A±(^ = 0) = 7°, (6.5c) 

where 7° = 7° = 7/, 7° = 7° = 7^. 

There is a good intuitive reason to separate the equa- 
tions for the coupling constants for i = 1, 4 from those for 
i = 2, 3. The latter group contains the quadratic interac- 
tion breaking the supersymmetry. Moreover, this would 
be the only group if we did not introduce the hermitiza- 
tion procedure of Sec. IV A. Those are the modes that 
will directly contribute to the observable quantities, see 
next Section, and we will call this sector "physical" . 

The coupling constants with z = 1,4 arc related to 
the fields that appear as a result of the Hermitization 
procedure. Their quadratic parts remain supersymmetric 
and that is why this sector by itself does not contribute 
to any observables. This is the reason why we will call 
this sector "non-physical". 

In the one loop approximation, Eq. (6.5) these two 
sectors do not talk to each other and we will consider 
them separately. 



C. Solution of RG equations in the "physical" 
sector. 

Equations (6.5b) is the system of the first order non- 
linear equations with the triangular structure (there is 
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no feedback of constants A, /? into evolution of the four 
particle vertex 7). As such, it can be easily solved with 
the initial conditions, Eqs. (6.5c): 



73(0 



(0 = 



(0 



C . 



(6.6) 



A+- (0 = A3-+ (0 



4* 



where we introduced the notation 



> 0. 



7.° 



(6.7) 



is defined in 



and the backscattering amplitude 
Eq. (4.49). 

From Eq. (6.6) we see that the forward scattering am- 
plitude is not renormalized in contrast to the backscat- 
tering ones. As wc consider the repulsive case, the am- 
plitudes 73, /Ja, 63 tend to zero in the limit T,\0\ ^0 and 
this is a "zero charge" situation. This behavior should, in 
principle, be seen using the conventional diagrammatic 
analysis. However, as the spin excitations considered 
here correspond to two particle electron Green functions, 
one should consider four particle electron Green functions 
in order to identify the behavior described by Eq. (6.6). 
Such interactions in four particles Green functions have 
not been studied previously and the result of Eq. (6.6) 
is a major and decisive step in this direction. The "zero 
charge" behavior does not indicate any drastic changes 
of the ground state of the system but is interesting on its 
own because it is definitely not present in the orthodox 
Fermi liquid picture. 

As the result corresponds to the zero-charge flow, the 
one loop renormalization group and Eq. (6.6) would solve 
the problem completely. The renormalized amplitudes A 
can be used for calculating the thermodynamic properties 
of the system, as it will be done in the next section. On 
this route, however, a potential reef may rise and we turn 
to the statement of this problem now. 



D. Solution of RG in "non-physical" sector and 
possible instability. 

Solving Eq. (6.5a) with the initial conditions (6.5c) 



71 (0 = Pt (0 = 7/; 



weobtain 



71 (0 = 
01 (0 



{6. 



where we introduced the notation 



7/(^) 



>o, 



(6.9) 



and forward the scattering amplitude 7° is defined in 
Eq. (4.49). 

The behavior of the amplitude 71 (^) from Eq. (6.8) 
comes as a real surprise because it demonstrates the ex- 
istence of a logarithmic pole. This pole should be reached 
at ^ = ^^'^ signal on an instability of the ground 
state because, at first glance, the scenario looks similar 
to the one leading to the BCS theory of superconduc- 
tivity. Does the logarithmic pole in Eq. (6.9) lead to a 
phase transition? We do not try to answer this question 
in this paper but the situation looks more complicated 
than usual because, within the RG scheme, the scatter- 
ing amplitude 71 does not enter thermodynamic quanti- 
ties like e.g. the specific heat, see the next section, and 
it does not couple to the physical sector at least on the 
level of one loop renormalization group. 

We can envision two different scenarios that may follow 
from the existence of the pole in 71 in Eq. (6.9): 

1. The amplitude 71 (^) does not enter any physical 
quantity and, therefore, the pole does not mean 
anything. In this case, we would be able to use 
all the equations for the backscattering amplitudes 
down to ^ = 0, which would allow us to go in tem- 
perature down to T = 0. This would mean that 
there are non-analytical corrections to the Fermi 
liquid but otherwise the Landau theory of Fermi 
liquid is a correct low-tcmpcrature limit. 

2. The logarithmic pole means that at ^ = a re- 
construction of the ground state (either in a form 
of phase transition or sharp crossover) occurs at a 
critical temperature T^, which would manifest it- 
self in a formation of the gap in the non-physical 
sector, that would affect physical degrees of free- 
dom. In this case the ground state would change 
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and one would expect formation of an order pa- 
rameter. Then, the present RG treatment of the 
effective action would be applicable for T > Tc 
only. The formation of the gap means a break- 
down of the Fermi liquid description for the spin 
excitations. As concerns the charge degrees of free- 
dom, we have seen that they decouple from the spin 
ones on the very early stage, see Sec. Ill, and the 
Landau kinetic equation is applicable for them for 
arbitrary small temperatures. 

It is still to be seen which of these two scenarios corre- 
sponds to the original action and we relegate this ques- 
tion to the further study''^. We emphasize, however, that 
the answer may differ for d = 1 and d > 1, because there 
is a large set of diagrams logarithmic in ID and non- 
logarithmic for d = 2.3, see e.g. Fig. 10. 



the contribution of the Cooper channel in VII D. For one 
dimensional systems approximations of Sees. VII A, will 
turn out not to be sufficient and we will have to take ad- 
ditional terms into account specific for one-dimensional 
systems, see Sec. VII D. 



A. General formulae. 

Using the diagrammatic method of the calculations one 
can always cut one of the Green functions and express the 
thermodynamic potential Q (T) in terms of a sum over 
bosonic Matsubara frequencies w„ = 27rTn transmitted 
through this particular Green function: 



(T) = T^MK) 



(7.3) 



VII. SPECIFIC HEAT 

We have performed the renormalization group calcu- 
lations for the case when the vectors n and n' of two 
spin excitations were close to each other (parallel or an- 
tiparallel motion). Only in this limit one obtains large 
logarithms that determined the renormalization of the 
vertices. A crucial question is whether or not this narrow 
region of the phase space can bring an important contri- 
bution to thermodynamics or other physical quantities. 
This is not quite trivial question because the system was 
not assumed to be one- or quasi-one-dimcnsional, and 
one could imagine that all the effect of the singularities 
in the vertices would be washed out after the summation 
over the whole phase space. 

In fact, this almost parallel motion of the spin excita- 
tions does not contribute much into the thermodynamic 
potential fig (T) itself. Fortunately, this is not a very in- 
teresting quantity and what one would like to know are 
derivatives of the thermodynamic potential with respect 
to temperature or other sources. In the present paper, 
we restrict ourselves with the specific heat 



C 



(7.1) 



Our goal is to identify the non-analytic contributions to 
the specific heat using the properties of the effective the- 
ory established in the previous sections. 

As we will show, our low-energy field theory is appli- 
cable for the calculation of 



Sfls (T) = (T) - (T = 0) , 



(7.2) 



and we will focus in this section on the calculation of the 
latter quantity. 

We will present the main approximations and manip- 
ulations suitable for any dimensions in Sec. VII A. We 
will collect the final results for two- and three- dimen- 
sional systems in Sec. VII B and discuss their relation to 



where M (w„) is a function of the frequency to be calcu- 
lated later. 

The sums of the type (7.3) are very often divergent at 
high frequencies if one uses expressions available from a 
low energy effective theory. This problem can be avoided 
calculating the quantity (T) from Eq. (7.2). Using the 
Poisson formula, we represent 5Cl (T) in the form 



(T) = ^ / M(a;)exp(- 



ilu) 



djjj 
2^' 



(7.4) 



which improves the convergence significantly. The essen- 
tial frequencies in Eq. (7.4)) are of the order of T and 
are smaller then those frequencies that form logarithms 
in the vertices. That is why the renormalized vertices 
calculated in the previous section become useful. 

To proceed with actual calculation, we notice that if 
we kept in the action S {-0} only the supersymmetric 
part Sa {V--} +^4 {V'}, see Eqs. (4.44), (4.45), (4.55) only, 
we would obtain unity for the partition function and, 
thus, no contribution to the thermodynamic potential O. 
The interaction terms S2 {V'} 'Uid S-^ {V^}, see Eqs. (4.51), 
(4.59), (4.52), (4.62) violate the supersymmetry and, as 
a result, one obtains finite contribution to Og only when 
expanding in such terms. 

As all the high-frequency lo > T contributions are al- 
ready included into the renormalized value of the vertices, 
the thermodynamic potential f2 (T) can be expanded in 
terms of the renormalized action S2 {^'}, Eq. (4.62), and 
the lowest non-vanishing orders take the form 



Qs (T) = Oi (T) + ^2 (T) ; 
Oi (T) = T(52{V})o; 

il2(T) = -|([52Wf)^: 



(7.5) 



and (...)q was defined in Eq. (5.2). The corresponding 
diagrams are depicted in Fig. lla,b). 

It will turn out that correction f^i is analytic function 
of temperature whereas the most interesting term, fl2 is 
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d) 




f) 




= 



FIG. 11: Lowest order corrections to the temperature depen- 
dent part of the thermodynamic potential SQ{T). The nota- 
tions are explained in Sec. V, see Figs. 4 - 9. The coefficients 
in the brackets correspond to the number of symmetries of 
the diagram. The dotted line for the Green function means 
that summation over the frequency of this line is performed 
with Eq. (7.4) replacing the summation (7.3). Diagrams c,d) 
vanish due to the supcrsyrnrrictry. Diagrams e,f) are sep- 
arately of higher order in the forward scattering amplitude 
and, moreover, cancel each other. 



non-analytic. For the pedagogical reasons, we will con- 
sider an analytic correction first and then use the gained 
knowledge for more involved analysis of the non-analytic 
contribution. 



1. Analytic contribution. 

Using the rules of the diagrammatic technique of Sec. V 
we can obtain analytic expression for fli (T) (for unit vol- 
ume) in terms of the function K introduced in Eq. (5.52): 



X ^ =(a;,k,n;u)A^i'"=(6' = 0;M,u;k); 

<Tl<72=± 

kj^ = k — n(k • n), rj^ = r — n(r • n) 



Ar"=(^;M,tt;k^) = /d'^re-^'^"-Ar'^Me;w,«;r^)/(|r|), 

(7.6) 

where the cutoff fimction / (r) is defined by Eq. (2.29), 
all the other entries were introduced in Sees. IV, V, and 



the convention (4.13) for the angular integration is used. 
Replacing the summation over the Matsubara frequency 
with the integrations according to Eqs. (7.3)-(7.4), and 
using Xij from Eq. (4.53), we obtain 



sni (T) 



3 / udu 
'0 



dn 



00 

■E 

=1 



duj 

2^^"P 



■ilui 



/(k) dk 
(27r)'' 



wpkn 



tii^kn 



A+- (^ = 0), 



(7.7) 



and the value of A^ is given by Eq. (6.6). The integral 
in Eq. (7.7) contains integration over all directions of 
the unit vector n with the normalization of Eq. (1.2). 
As we obtained in the previous section, Ag ~ is not re- 
normalized by interaction and it is given by its bare value 
7/ = ^j{9) = 0. Nevertheless, the remaining integral in 
Eq. (7.7) gives a temperature dependent contribution and 
let us show how one can calculate this integral. 

The integration over u is trivial and gives 1/2. The 
integration over k can be performed separately for the 
component k\\ parallel to n and kj^ perpendicular to n. 
Essential fcy are of the order of T/vp, whereas k±_ are of 
the order of the maximum momentum kc — . There- 
fore, we can neglect fc|| in the cutoff function. Integration 
over k]\ is, then, immediately performed with the result 



5VLi (T) 



lim 

??^-l-o 

37/ 
2vf 



1=1 ■' 



dio 



.lu) 



— |a;|expi -I— -r/iwl 



/o {k±_ro) 



d<^- 



(27r) 



d-l 



The small parameter ry in the exponential is added to 
provide the convergence of the integral over oj. 

Expression in the first line of Eq. (7.8) takes after the 
w-integration the following form: 



1 



lim 5- + 



{■KTf 



(7.9) 



Substituting Eqs. (7.8)-(7.9) into Eq. (7.1) we obtain the 
corresponding contribution 5C\ (T) to the specific heat 



5Ci 



^(37/) 



/o {k±ro) 



fi'^-ik_L 
(2^ 



d-l ■ 



(7.10) 



The parameter Aq — 1/ro becomes of the order of kp^ on 
the limit of the applicability of the theory. 

We see from Eq. (7.10) that the correction SCi (T) does 
not change the linear dependence of the specific heat on 
temperature. 

Let us discuss the significance of this result and its rela- 
tion to the free mode consideration of Sec. Ill C. First of 
all, direct comparison shows that Eq. (7.7) is equivalent 
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to the first term in the pcrturbativc expansion of temper- 
ature dependent parts of Eqs. (3.34'), (3.30'). Therefore, 
the singlet contribution (3.30') also has to be taken into 
account in Eq. (7.10), which leads to the replacement 



(37/) ^ 37/ 



where 



7/^ 



Y{o = 0). 



Equation (7.10) gives the contribution which does not 
depend on the cutoff in one dimension. In this case it 
describes the rcnormalization of the velocities of the spin 
and charge modes in Luttinger liquid regime. In con- 
trast, in higher dimensions, d > 1, the coefficient does 
depend on the cut-off. However, all this contribution 
coming from the small distances can be ascribed to the 
rcnormalization of the effective mass and included into 
the partition function of non-interacting quasiparticles, 
[Zq from Eq. (3.12)]. 

The effective theory of the interacting collective modes 
being the effective low energy theory does not describe 
such ultraviolet corrections and that is why we cannot 
identify the numerical coefficient from our theory. How- 
ever, the effective low energy theory'^" does describe the 
non-analytic corrections as the latter are associated with 
the spatial scales ~ vp/T l/kp- We turn to the 
derivation of such non-analytic corrections now. 



2. N on- analytic contribution. 



The analytic expression for the diagram Fig. lib reads 
[see Eq. (7.6) for the notations] 



oj„ i=2,3j=l 



du\du2 / / dnidn.2 



" ^ ' CTl ,...,Cr4=± 

X {kJ^'^^ (w, k, m; ui)Af^<^^ (nTS^; ui, «2; kj 
X Kf^* (w, k, n2 ; W2) A^^-^i {n^i ; uj, «i ; k^) } ; 



(7.11) 



where the convention (4.13) for the angular integration 
is implied, and in definition of k^, the direction n means 
n = (ni + n2) /2. Using Eqs. (5.52) and (4.53), keeping 
only terms with the pole location such that the result 
does not vanish after integration over k» = k • n, and 



'^'^ Without any drawbacks one can formulate the effective low- 
energy theory for the Fermi liquid function such that 87/ = 7^- 
The ultraviolet correction to the specific heat does not appear in 
this case at all. 



using replacements Eqs. (7.3)-(7.4), we find 

diu / ,lu} 



6^2 (T) = - 
X Idnidn2 



6 lim > 



(2^^"PI"'t 



d'^k 



Y (rrjlii;k_L; fcy) (u;,k; ni,n2) , 

(7.12) 



and the convention (4.13) for the angular integration is 
implied. The function Y {6; \s.±) defined as 

Y [0;'k±,k\\) = jj uiU2duidu2 
^ { [^3"" (6';wi,M2;k_L,fc||)]^ 

+ (6*; Ml, W2; k_L, fc||) Ag'" (6*; ui, U2; k_L, | 

(7.13) 

will be the most important entry in the final expression 
for the specific heat. In Eq. (7.13), we wrote explicitly 
the transverse and the longitudinal momenta of k, see 
Eq. (7.7), as their role will be different. 
The formfactor 



?rf (w,k;ni,n2) 



{iuj + ?;i?kn2) {iuj — ypkni) 
{iuj — VFkn2 ) {ico + vpkni) 



(7.14) 



depends on the dimensionality of the system and it de- 
scribes basically the free propagation of the two spin ex- 
citations in almost opposite directions. 

If one used the bare values of A, see Eq. (4.65), we 
would obtain the second term in the expansion of Eq. 
(3.34') in powers of 7. The main advantage of Eq. (7.14) 
is that it accounts for the logarithmic rcnormalization 
of the quadratic interaction obtained in Sees. V,VI, see 
Eq. (6.7). 

The non-analytic contributions originate from the 
small region of the phase space |ni — n2| <C 1. It en- 
ables us to introduce [cf. Sec. V B] 



n = (ni -I- n2)/2; 5n = ni - n2; 
fell = k • n; k_L = k — fciin, 



(7.15) 



and integrate over fcy in Eq. (7.13). 

To facilitate the integration we introduce the function 

f diodkn f Ilo \ 
(5n; fcx) = ^hm^ ^ J ^ exp (^-^- -rj\u;\j 



^ 2^ J'd ( ±kj_ -I- fcyn; n -I- — , n - 



(7.16) 



Using 



00 « 
lim V / 



cL;exp ( -i^ -?7|w| 



J dfc||/2(k)^ =0, 
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we obtain 



Vd=i = 0, 



(7.17a) 



which means that for one dimensional systems there is 
no contribution to the specific heat of the second order 
in the backscattering ampHtude in accord with known 
results^^. This does not mean, however that there arc no 
logarithmic contributions to the specific heat in Id at all. 
Corrections of the third order in the effective amplitude 
will be recovered in Sec. VII D. 

For the higher dimensions we obtain 



T f dcv f Ilo 

^<^=2.3 = ^lnn^l^y ^exp(^-^-. 



7?|W 



X W - 



4a;2 + (5k_Ln)^ 

(7.17b) 

The first term in Eq. (7.17b) is similar to Eq. (7.8) and 
produces only analytic contribution to the specific heat. 
We will disregard this term and focus only on the second 
contribution. This is this contribution that quantifies the 
effect of the small region of the phase space where the 
interaction renormalization is strong. It is worth not- 
ing that the characteristic k\\ contributing into the result 
where of the order of T /vp and, therefore, the separa- 
tion of the integration on the transverse and longitudinal 
parts is well justified. 

Substituting Eq. (7.17b) in Eq. (7.13) and performing 
the Fourier transform in the transverse direction we find 
for d > 1 



hm y / 



2tt 



■ cxp 



2 \ujx\ iiol 



(7.18) 



where the integration over n and 5n are performed using 
conventions (4.13) and (1.2) respectively, and 



^(^;|r^l 



d'^-^k 



£f^e*--y(e;k^;k||=0). 



(2^) 



(7.19) 

Integration over co and n in Eq. (7.18) can be immedi- 
ately performed with the result 



xY{\Sn\; x)Re | ^.rp. 



(7.20) 



where integration over Sn has to be understood as inte- 
gration over usual d—1 dimensional vector with \dn\ < 1. 



The remaining integrations in Eq. (7.18) arc slightly 
different for c? = 2, 3 and we will describe them sepa- 
rately. For d — 2, 6n is a. one-dimensional variable. After 
obvious rescaling, Eq. (7.20) gives 



snf (T) 



OO 



dx 



X Re 



I 



# 41^ Y 



(7.21) 



The integral in Eq. (7.21) convergent at ~ 1. Because 
a; ~ ro, and T <^ vp/ro, the upper limit in the integral 
can be put to infinity. On the other hand, according to 
Eq. (6.4), Y{0) does not depend on 6» at 6» < Tro/vp, 
therefore, we can put the first argument in Y equal to 
zero. The remaining integral takes the form 



f 

Jo 



d<i)(lP 



2' 



/o {icl> + ly 
The sum over I is trivial, and we obtain 



6nr (T) 



6C (3) r3 



TTVt 



Y{9 = 0)., 



(7.22) 



where C(3) « 1.202 ... is the Riemann ^-function and 

Y {6) = Y (9;^ = 0) , (7.23) 

see Eq. (7.13). The value 6* = in the function Y {9) 
entering Eq. (7.22) corresponds to the exactly backward 
scattering. 

Let us consider, now the three-dimensional case. The 
result turns to be logarithmic with the main contribution 
from Tx/vp ^ \Sn\ ^ 1- Expanding the last factor in Eq. 
7.20and keeping only the first non-vanishing contribution 
we find 



(T) 



36T^ ^ 1 



(fdn 
|(5n|2 



xdxY {\6n\; x), 
(7.24) 



As the integrand decays rapidly at a; ~ ro, the upper 
limit in the last integral can be put to infinity for 6 = 
\5n\ > Tro/vp. Using 



27r jxdxYiO; x) 





Y{9), 



see Eq. (7.23), we obtain 



sn't^ (T) 



27r2T4 r 
'W J 

Trp 



d9 Y {9) 



(7.25) 
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Equations (7.22) and (7.25) arc the main results of this 
subsection. They show that the thermodynamic poten- 
tial for the both two and three dimensional systems can 
be expressed in terms of the function Y (6) , that has to 
be determined from Eqs. (7.13), and (6.6). Carrying on 
this program and obtaining the final expressions for the 
specific heat will be subject of the next subsection. Here, 
we just emphasize the difference between Eqs. (7.22) and 
(7.25). Namely, in 2d case one has to take the function 
Y at strictly 9 = 0, whereas in 3d it involves integrals 
over all angles. At the same time, even in 3D, the main 
contribution comes with the logarithmic accuracy from 
small angles 9 < 1 , that describe again the scattering 
close to backward one. 



B. Final results for two- and three- dimensional 
systems. 

From now on we will concentrate only on the non- 
analytic contribution Sil2~^'^ that is why we will 
omit the subscript. Before writing the final results in 
a general form, it is instructive to perform the calcula- 
tion replacing the quadratic interaction constants A^^ 
in Eqs. (7.13) and (7.23) with their bare values (4.65). 
After integration over ui,2 in Eq. (7.13) one obtains 

Y^'H0) = lbbi9)ffi0) = llU0), (7.26) 

where ■yb{9) is defined in Eq. (4.49), and / (k) is given 
by Eq. (2.13). Substitution of Eq. (7.26) into Eqs. (7.22) 
and (7.25) yields 



C(3)T3 



TTVz 



76 = 76(^ = 0), 



(7.27a) 

where we restored the explicit subscript "t" for the 
triplet contribution. If Y^^^ (9) is a smooth function for 

9^0, the three-dimensional thermodynamic potential 
reads with the logarithmic accuracy 



T 



(7.27b) 



Similarly to Eq. (7.10), formula (7.27) is equivalent to 
the second term in the perturbative expansion of temper- 
ature dependent parts of Eqs. (3.34'), (3.30'). Therefore, 
the singlet contribution (3.30') has also to be taken into 
account in Eq. (7.27), which leads to the replacement 

(376^) ^ (3^1 + where Yb ^ Y{6 = ^r). Unlike 

the analytic contribution, Eqs. (7.27) do not contain the 
ultraviolet cut-off and their contribution can not be as- 
cribed to the renormalization of the effective mass. These 
are clear effects of the contribution of the bosonic collec- 
tive modes of the Fermi liquid into the thermodynamics 
of the system. 



Using Eqs. (7.27) in Eq. (7.1) and introducing the den- 
sity of particles 

we write the correction to the specific heat per particle 
5c = 5C/N as 



6c 



(bare) _ 3C (3) f T \ 



d=2 



-) {H + btf)-, (7.29a) 



5c 



(bare) 
d=3 



(7.29b) 

[The last term in Eq. (7.29) is the contribution of the 

third order in coupling constant which was obtained in 
Ref. 28 and we refer the reader to this reference for the 
explicit form of the c;ocfficients in this term and will not 
write explicitly in the subsequent considerations.] It is 
worth reminding that pi? is not renormalized by interac- 
tion, whereas tp ^ vppp/'^ is significantly affected^. Ac- 
tually, vp has the meaning only as a quantity describing 
the slope in a leading linear in temperature quasiparticle 
contribution to the specific heat. 

Equations (7.29) agree with the corresponding expres- 
sions obtained previously in a number of works using con- 
ventional diagrammatic expansions (see e.g. Refs. 18,28, 
for the latest developments consisting in accurate evalu- 
ation of the angular and q integrals in expressions similar 
to Eqs. (3.30), (3.34) and obtaining correct analytic ex- 
pressions for the first time) '^™. 

Using the conventional diagrammatic technique one 
can hardly go beyond the first orders, which would be def- 
initely enough for the singlet channel. At the same time, 
using the present bosonization scheme we have found for 
the first time the logarithmic contributions discussed in 
the previous Sections and have derived and solved proper 
renormalization group equations. Now we can include the 
logarithmic contributions into the formulas for the spe- 
cific heat. The only thing that remains to be done is to 
calculate the function Y {9) from Eqs. (7.23), (7.13). 
Using the explicit expressions (6.6) for the couplings 

3 



A^^, and the formula 



/i" 



duidu2 {U1U2Y 



(1 + xiiiiiL-zY + (1 + -r-iiLiUaY 
{I + X1U1U2) {I + X2U1U2) 



{Xi + X2)uiU2 
2{1 + XlUlU2)^ (1 + X2UlU2)^ 
1 



2(1 + Xi)(l+X2)' 



The definitions of the coupUng constants 7, Eq. (3.35), diff'crs 
by the factor 2 from the one of Ref. 18,28 because we used the 
density of states v per one spin orientation but not the total 
density of states. 
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we reduce Eq. (7.13) to the form 



(^) ^w. 



(7.30) 



where function /j_ {r±/ro) is given by Eq. (5.16b). The 
variable X is defined as 



Xie) = -pidlM In [max {9, T/e^}] ; 
X(r) = X(^ = 0) = M<i75ln(^), 



(7.31) 



with the parameter fXd given by Eq. (5.16), and the cut- 
off energy defined as, Eq = ^ ^ e p ■ 

Equation (7.30) gives the most general form of the 
function Y {9) for any function /. The asymptotics of 
the function Y {9) in the limit X -C 1 can easily be writ- 
ten as 

Y {0) = m - mrt' I f (k^) . 

(7.32) 

where the cutoff function / (k) is given Eq. (2.13). 

The first term in Eq. (7.32) is what has been used 
when deriving Eqs. (7.29). The second term leads to an 
additional logarithm in the T'^-dependence of the spe- 
cific heat in both two and three dimensional systems. 
Notice, that the latter term depends on the form of the 
ultraviolet cut-off /(k). As the function / (k) has been 
introduced in our theory phcnomcnologically, the com- 
plete results depend on its form, which is a deficiency 
of our low energy bosonization approach, where the bare 
coupling constants were introduced independently of the 
ultraviolet cutoff. At the same time, we will sec that it 
gives at least a good qualitative description of the inter- 
esting temperature behavior, and will address the issue of 
the dependence of the coupling constants vs. the cutoff 
function in the next subsection. 

A simple expression for the specific heat for arbitrary 
X can be obtained if we choose the function f± {r±/ro), 
Eq. (5.16a) as 

h {r±/ro) = exp (-r^/ro) (7.33) 
The choice (7.33) in the coordinate space corresponds 



to 



(l + k2r2) 



2^3/2 



(7.34) 



in the momentum space. 

Substituting Eq. (7.33) into Eq. (7.30), and perform- 
ing the remaining integration, we obtain 



Yd=3 {9) 



(7.35) 



2[X{e)] 



2 ' 



where Li2{x) = 'Yl'kLi^^ 1^"^ ^^e polylogarithm fimc- 
tion. Using Eq. (7.35) we can write the asymptotics of 
the function Y (6») for X > 1 as 



ld=2 {9) 



2X2 ' 



lb 



2ln*X 



8X2 



(7.36) 



The asymptotic behavior X^^ {^i Eq. (7.36) is not very 
sensitive to the form of the function / (k), although the 
power of In X is not universal. On the other hand, more 
accurate treatment of those double logarithmic temper- 
ature dependencies would be an overstepping of the ac- 
curacy of the one loop approximation anyway. 

The final formula for the specific heat 5c per particle 
can be written with the logarithmic accuracy as 

3C(3)r2 [^^,^2 , 37g{ln[l + X(T)]}2 | 

(7.37a) 
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(7.37b) 



where the variable X(r) is defined in Eq. (7.31). Equa- 
tions (7.37) refine Eqs. (7.29) by including all the leading 
logarithmic corrections originating from the interaction 
of the spin modes. 

Equations (7.31), (7.32), (7.35), (7.36)) give the final 
results for different cases and demonstrate non-trivial log- 
arithmic dependencies on temperature. This behavior is 
more complicated than what one usually expects for the 
Fermi liquid picture. The unusual behavior is due to the 
interaction between the spin excitations. As concerns the 
charge excitations, they contribute in a more simple way, 
and their contribution is completely expressible in terms 
of the Fermi-liquid interaction function. 



C. On the role of the Cooper channel and choice of 
ultraviolet cut-off function. 

Strictly speaking, all the results we present here can 

be justified for tqPf ^ 1, where vq has been introduced 
as the shortest length of our low energy theory. How- 
ever, we hope that they remain relevant for the initial 
model of the Fermi gas with a repulsion. The scale ro in 
this case has to be found from the explicit calculation of 
the logarithmic corrections in original model of the inter- 
acting fermions rather than in the reduced model (2.10). 
Such calculation will not be done in the present paper, 
however, we will try to outline the steps which should be 
performed, without claiming too much rigor. 

It is well known^^ that the Fermi liquid functions expe- 
rience the strong logarithmic renormalization for scatter- 
ing directions close to backwards. (Such a logarithmically 
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divergent term for the Fermi liquid function has been ob- 
tained for the first time in Rcf. 59.) This is because the 
Cooper channel, see Fig. 12, at such angles has a strong 
mixing with the electron-hole channel. If the logarithmic 
renormalization were present only in the Cooper channel 
the result would read 



£ = ln 



70 



max (|w|, ufIQI 



(7.38) 



i.e, one had a zero-charge situation provided operator 70 
is positive definite. The needed function jc{9;Q,ui) are 
defined as kernels of the operator jc- Notice, however, 
that the result dependends on all of the angular harmon- 
ics of the bare interaction 70 . 

There is a region in the phase space, however, where 
the Cooper channel and the triplet channel cannot be 
distinguished from each other, see also Figs. 2, 3, and this 
is the region that was studied in previous sections. Even 
though the structure of the results (7.37), indicate that, 
indeed, the result can be factorized to the logarithmic 
renormailization of the backscattering amplitudes as it 
was argued in Ref. 18, however, it does not mean that 
they coincide with the renormalization by the Cooper 
channel only. 



D. Peculiarities for one dimensional systems. 



The purpose of this subsection is to find the lead- 
ing singular correction to the specific heat the in one- 
dimensional case and compare the result with the one 
obtained in Ref. 55 for the spin chain. The low-energy 
properties of the latter model are the same as for the spin 
dynamics for the interacting electrons. 

Wc represent the temperature dependent part of the 
desired correction as [cf. Eq. (7.4)] 



E 

1,(2/0 



. 27r 

X [Ra (^02) 
duJldL02 



h (W2) + Kc (W2)] 

exp I h 



(7.39) 



where the subscripts a, 6, c, e indicate the analytic expres- 
sion for the corresponding diagrams on Fig. 13. 



pi + |,ei + fi 



a"2 



Pi - i,ei 



P2 — ^ , £2 



Ic ; Pi -I- P2; ei -I- £2 

Pf 



P2 -h f , £2 



Pi Pi 




-Pl -pi 



P2 



-P2 



FIG. 12: Loading logarithmic renormalization of the vertex 
76. The integration over the intermediate momenta Pi has 
to exclude the region |pi — pi,2| < 1/ro as the latter has 
already been included in the effective energy theory for the 
spin excitations, see also Figs. 2, 3. 



Summation over the matrix indices is performed using 
the formula analogous to Eq. (5.52) with the result 



£a,6 = -- E ^"^"'K,m(^2); 



0"1.2=± 

K-(-.) = 2 1ini^ lim /f^/^ 

rji— >0 r)2^5)i J Ztt J [ Iiry 

X cos(fcir?i +A:2?72)3?'^"^'(t^i,a;2;A;i,2), 
7er-(a;2)= lim lim 

X cos(fci77i + k2rj2) ^^^^'^ {u>2 - wi, wi; fci,2); 

K^^H^2)= lim lim / 

X cos (fci?7i -I- ^2772) 'Ji"^"^ {uj2 - a;i,a;i; fci,2); 

(7.40) 



where Y"^"^ depends on the interaction constants only: 



Closing this section, we notice that the main contribu- 
tion in the Kohn-Luttinger^° scenario of the supercon- 
ducting instability also originates from the region of the 
phase space studied in our paper. As in this region the 
Cooper channel intervenes the particle-hole channel, the 
simple use of the second order screened interaction in the 
Cooper channel^" does not appear to be justified. 



F'^i'^^ = jj^ {uiU2f duidu2[l^l'"'^ lis"' fis"^] . (7.41) 



The most interesting factors ^R'^^'"^ (wi, W2; A;i,2) are given 
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UJl +Ul2, 

ki + k2 




UJ2 — U)l, ki 



'■ UJ2,k2 + ki 





d) 



e) 







'1 




= 








2 



a»2 




FIG. 13: Diagrams of the third order, 5Q3{T) oc 75 , giving the 
leading logarithmic contribution the in one-dimensional case. 
The notation are introduced in Sec. V, see Figs. 4-9. The 

cocfScicnts in brackets corresponds to the number of symme- 
tries of the diagram. The dotted hrie for the Green function 
means that the summation over the frequency of this line is 
performed with Eq. (7.4) replacing the summation (7.3). The 
solid line implies the simple integration over frequency as at 
T = 0. Diagram d) does not contribute to the specific heat. 



by 

Ol++ = 

01 
01 



{iuj2 + ^2) (iwi - fci) 

2^(wi,2; fcl,2) 
(iw2 + ^2) [i (wi + W2) - (^2 + ^1)] . 
I'('^i,2; A;i,2) 
+ _ [i {u>i + UJ2) + (fc2 + fci)] {ioji - fci) ^ 

2?('^i.2; fci,2) 

[i (Wi + W2) + (/C2 + fci)] (Wi + W2) 



(fc2 



(7.42a) 
(7.42b) 

(7.42c) 



25(^1,2; fci, 2) 



(7.42d) 



©(^1,2; fcl,2) = (iWl + Vpkl) {iLU2 - VFk2) 

X [i{uJi + LU2) - vpiki + fc2)] 
X [i(wi + UJ2) + vpiki + fc2)] . 

Finally, 

Mc = — // du\du2 (wiM2)^ 73Ar"'"Aj"~ lim lim 
V J Jq m^om^-ni 

I f dkidk2 f dw\ cos (fci?7i — fc2?72) 
J J {2-Kf J 27r [iuji - vpki) {iw2 + VFk2) ' 

(7.43) 

and all the other contributions Mc either vanish or pro- 
duce contributions independent on 0^2. 

As usual for one-dimensional system with the lin- 
earized spectrum, the integrals (7.40), with the inte- 
grands (7.42), and Eq. (7.43) have the anomalous char- 
acter: each term could be eliminated by the shifts of 
the momentum if such arbitrary shifts were allowed. 
That is why the chosen order of limits is very cru- 
cial for the complete definition of the action. In 
terms of the original model, it corresponds to the reg- 
ularization of the singular terms (0j^(a;)0i(a;))0jj(x), 
and {(l)\{x)(j)]^{x)){(j)^j^{x)<pfi(x)), etc. in notation of 
Eq. (4.10), by shifting the coordinate of the left (right) 
movers by r?i(— ?72)- Such shifts eliminate all the diver- 
gent terms. 

Calculating the integrals (7.40) we find 



1^2 1 . 

\C02\ 
'4^' 



n: 



7^+- 



sgna;isgnw2 



4^1 



■ 7^-+ 



f 7^ 



(7.44a) 

(7.44b) 
(7.44c) 
(7.44d) 



and 



l^2| 



'-^11 



du2{uiU2) isAg+A^^ . 

(7.44e) 



Substituting Eqs. (7.44) into Eq. (7.39), using 
Eq. (7.9), the explicit form of the coupling constants, 
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Eq. (6.6), Hi = 2, u = 1/{ttvf), and the formula 



ii: 



duidu2 {U1U2) 



(1 + XU1U2) 
1 



(1 + XU1U2) 



3(1 



we find the leading logarithmical contribution to the ther- 
modynamic potential 



III 
16vf 



276 



l + 27bln^ 



(7.45) 



This correction agrees with the one previously obtained 
one for the spin chains (see Eq. (3.17) of Ref. 55) and we 
conclude that the supersymmetric low energy theory de- 
veloped in this paper reproduces all of the known physical 
results despite the fact that the intermediate degrees of 
freedom apparently differ from those for the conventional 
bosonization. 

One can also prove by the explicit calculation that 
there is no contribution to the specific heat proportional 
to 7/7^. Corresponding diagrams are shown on Fig. 14. 



VIII. DISCUSSION 

We have considered thermodynamics of an electron gas 
with a repulsion in arbitrary dimensions. This model be- 
longs to the class of systems that should manifest a Fermi 
liquid behavior, which implies that the temperature be- 
havior of thermodynamic quantities is similar to the one 
of an ideal Fermi gas. 

In order to investigate low lying excitations like spin 
or charge ones we developed a new method of bosoniza- 
tion that allows us to replace the initial electron model 
by a model for low lying excitations. Our approach 
is based on a method of quasiclassical Green functions 
and differs from earlier high dimensional bosonization 
schemes^'^'^^"'*^ . In contrast to the latter approaches we 
can consider not only charge excitations but also spin 
ones. This advantage is crucial because the spin exci- 
tations arc much more interesting than the charge ones. 
The importance of the spin excitations is seen from Eqs. 
(3.10a, 3.10b). In contrast to the charge excitations, the 
spin ones interact with each other via the fluctuational 
magnetic field h. 

Studying the low lying spin excitations we have dis- 
covered non-trivial logarithmic contributions to ther- 
modynamic quantities originating from their interaction 
and succeeded in summing them using a renormalization 
group scheme. The logarithmic contributions come from 
momenta of the two excitations parallel or antiparallel 
to each other (forward [which docsnot contribute to the 
physical quantities studied in this paper] and backward 
scattering). To some extent, the system manifests one 
dimensional properties even if we work in two or three 
dimensions. 



a) 




d) 




FIG. 14: Diagrams of the third order, S^l'^-^' (T) oc j^^f. Those 
diagrams do not contribute to the specific heat. Non-filled 
vertices correspond to the physical forward scattering ampli- 
tudes which are not renormalized. Diagram b) vanishes due 
to the supersymmetry. Diagram a) vanishes separately. 



In principle, wc coidd solve Eq. (3.10b) using a pertur- 
bation theory in the effective field h. Then, substituting 
such a solution into Eq. (3.11) we would be able to cal- 
culate the thermodynamic quantities. However, this is 
not a completely safe scheme. The problem is that the 
linear operator acting on the variable Sn is Eq. (3.10b) 
is not Hermitian because it contains linear derivatives in 
time and coordinates. It is well established^^'^^ that lin- 
ear derivatives in the operators can lead to a new physics 
because, the non-hermitian operators may have complex 
eigenvalues. This problem can be avoided by the process 
of the hermitization^^ and everything can be reformu- 
lated in terms of a field theory containing supervectors. 



41 



This is why wc developed a supersymmetric scheme and 
used it for the calculations. This gives also advantages 
because wc can demonstrate the renormalizability of the 
theory, which is difficult using the perturbation theory. 

The method we developed is applicable also for one 
dimensional systems and it exactly reproduced the known 
result for the logarithmic correction to the specific heat. 
It gave us a great deal of confidence in the correctness of 
our procedure, However, our main interest is the higher 
dimensional systems and we do not intend to compete 
with the very well developed methods in Id. 

Using the method of the renormalization group we 
calculated all relevant vertices of the theory, which al- 
lowed us to calculate the thermodynamic potential and 
the specific heat, Eqs. (7.37a, 7.37b). In the lowest order 
one neglects the interaction between the spin excitations 
and obtains Eqs. (7.29a, 7.29b) that have been obtained 
previously by conventional diagrammatic expansions (see 
the latest works^*'^^ and references therein). These cor- 
rections are already non-analytic in T^, which was the 
main motivation for their previous study. We see from 
the results obtained here that the problem is even more 
interesting and the temperature behavior of the thermo- 
dynamic quantities in really non-trivial. 

We derived the results in the weak coupling limit and 
the approximations we used are justified. Although we 
cannot apply the results in the strong coupling limit, it is 
difficult to imagine that the non-trivial temperature de- 
pendence would not be relevant in that region. This can 
lead to complicated effects near quantum critical points. 

The application of the RG scheme we developed has 
lead us to a very unusual result, namely, the amplitude 
7i describing the forward scattering of the spin excita- 



tions has a logarithmic pole, Eq. (6.8), and can diverge 
below a critical temperature. We emphasize that the 
appearance of this diverging vertex is a consequence of 
the hermitization procedure we used for the derivation of 
the field theory and its existence cannot be noticed using 
conventional diagrammatic expansions. 

Very often divergencies of scattering amplitudes lead 
to a phase transition as, e.g., in theory of superconduc- 
tivity. At the same time, the logarithmic pole in the 
Kondo problem does not lead to any phase transition. 
The situation now is even more tricky because the for- 
ward scattering amplitude does not enter the thermody- 
namic quantities in the perturbation theory and the RG 
scheme at all. Therefore, even if it diverges, this does not 
necessarily mean a phase transition because it may drop 
out from physical quantities. Clarifying this situation is 
the most challenging continuation of the present study. 



IX. ACKNOWLEDGMENTS 

We enjoyed interesting discussions with A. Chubukov, 
G. Schwiete and R. Teodorescu. We are grateful to F. 
Essler and A. A. Nersesyan for bringing Ref. 55 to our 
attention. We are grateful to D. Basko, G. Catelani, 
A, Chubukov, L. Glazman, D. Maslov, and A. Millis for 
reading the manuscript and valuable remarks. K.B. Efe- 
tov is thankful to Columbia Physics Department for the 
kind hospitality and to grant Transregio 12 "Symmetries 
and Universality in Mesoscopic Systems" of the German 
Research for financial support. 



1 L.D. Landau, Zh. Eksp. Teor. Fiz. 30, 1058 (1956) (Sov. 

Phys. JETP, 3, 920 (1956)); ibid 32, 59 (1957) (Sov Phys. 

JETP 5, 101, (1957)). 
^ T. Timusk and B. Statt, Rep. Prog. Phys. 62, 61 (1999). 
^ M. Norman and C. Pepin, Rep. Prog. Phys. 66, 1547 

(2003). 

* P.W. Anderson, Phys. Rev. Lett. 66, 3226 (1991). 
^ R. Shankar, Rev. Mod. Phys. 66, 129 (1994). 
® W. Metzner, C. Castellani and C. Di Castro, Adv. Phys. 
47, 317 (1998). 

A. Houghton, H.-J. Kwon, and B. Marston, Adv. Phys. 

49, 141 (2000). 
® J. Kim and D. Coffey, Phys. Rev. B 62, 4288 (2000). 
® A. Abanov, A. Chubukov, and J. Schmalian, Adv. Phys. 

52, 119 (2003). 

^° L.D. Landau, Zh. Eksp. Teor. Fiz. 35, 97 (1958) (Sov. 
Phys. JETP, 8, 70 (1959)). 

A. A. Abrikosov, L.P. Gorkov, and I.E. Dzyaloshinskii, 
Methods of Quantum Field Theory in Statistical Physics, 
(Prentice Hall, New York, 1963). 

D. Pines and P. Nozieres, The Theory of Quantum Liquids, 
(Benjamin, New York, 1996). 

A.L. Fetter and D. L. Walecka, Quantum Theory of Many- 



Particle Systems, (MacGraw-Hill, New York, 1971). 

G. Baym and C. Pethick, Landau Fermi liquid Theory, 

(Wiley, New York, 1991). 

G.M. Ehashberg, Sov. Phys. JETP 16, 780 (1963). 

S. Doniach and S. Engelsberg, Phys. Rev. 17, 750 (1966). 

W.F. Brinkman and S. Engelsberg, Phys. Rev. 169, 417 

(1968). 

A.V. Chubukov, D.I. Maslov, and A. Millis, cond- 

mat/0509163. 

D.J. Amit, J.W. Kane, and H. Wagner, Phys. Rev. 175, 
313 (1968); ibid 175, 326 (1968). 

D. Belitz, T.R. Kirkpatrick, and T. Vojta, Phys. Rev. B 
55, 9452 (1997). 

D. Coffey and K.S. Bedell, Phys. Rev. Lett. 71, 1043 

(1993). 

A.V. Chubukov and D.L. Maslov, Phys.Rev. B 68, 155113 

(2003) ; ibid 69, 121102 (2004). 

V.M. Galitski and S. Das Sarma, Phys. Rev. B 70, 035111 

(2004) . 

G. Catelani and I.L. Aleiner, Zh. Exp. Teor. Fiz. 127, 372 

(2005) ; [JETP, 100, 331 (2005)]. 

M.A. Baranov, M.Yu. Kagan, and M.S. Mar'enko, JETP 
Lett. 58, 709 (1993). 



42 



^'^ G.Y. Chitov and A.J. Millis, Phys. Rev. Lett. 86, 5337 

(2001); Phys. Rev. B64, 054414 (2001). 
^'^ A.V. Chubukov, C. Pepin, and J. Rech, Phys. Rev. Lett. 

92, 147003 (2004). 

A. V. Chubukov, D.L. Maslov, S. Gangadharaiah, and L.L 
Glazman, Phys. Rev. B 71, 205112 (2005). 

^® L.P. Gorkov, A.I. Larkin, and D.E. Khmelnitskii, Pis'ma 
Zh. Eksp. Tcor. Fiz. 30, 248 (1979) [Sov. Phys. JETP Lett. 

30, 228 (1979)]. 
^° E. Abrahams, P.W. Anderson, and T.V. Ramakrishnan, 
Phil. Mag. 42, 827 (1980). 
F. Wegner, Z. Phys. B 35, 207 (1979). 
K.B. Efetov, Adv. in Phys., 32, 53 (1983). 
K.B. Efetov, Supersyrnmetry in Disorder and Chaos, 
(Cambridge Univ. Press, New York, 1997). 

B. A. Muzykantskii and D. E. Khmelnitskii, Pis'ma Zh. 
Eksp. Theo. Fiz. 62 68 (1995) [JETP Lett. 62 76 (1995)]. 
K. B. Efetov, V. R. Kogan, Phys. Rev. B 67, 245312 
(2003). 

A.M. Finkel'stein, Zh. Eksp. Teor. Fiz. 84, 168 (1983) (Sov. 
Phys. JETP 57, 97 (1983)); A.M. Finkel'stein, in Electron 
Liquid in Disordered Conductors, edited by I.M. Khalat- 
nikov, Soviet Scientific Reviews, Vol. 14 (Harwood, Lon- 
don, 1990). 

^"^ G. Schwiete and K.B. Efetov, Phys. Rev. B 71, 134203 

(2005). 

3* A. Luther, Phys. Rev. B 19, 320 (1979). 

^® F.D.M. Haldane, Hclv. Phys. Acta. 65, 152 (1992). 

*° A. Houghton and J.B. Marston, Phys. Rev. B 48, 7790 

(1993); A. Houghton, H.-J. Kwon, and J.B. Marston, Phys. 

Rev. B50, 1351 (1994). 
*i A.H. Castro Neto and E. Fradkin, Phys. Rev. Lett. 72, 

1393 (1994); Phys. Rev. B 49, 10877 (1994). 

P. Kopietz and K. Schonhammer, Z. Phys. BlOO, 259 

(1996). 



P. Kopietz, Bosonization of Interacting Fermions in Arbi- 
trary Dimensions, ( Springer Verlag, Berlin, 1997). 
D.V. Khvcshchcnko, R. Hlubina, and T.M. Rice, Phys. 
Rev. B48, 10766 (1994). 

D.V. Khveshchenko, Phys. Rev. B 49, 16893 (1994); ibid. 
52, 4833 (1995). 

C. Castellani, C. Di Castro, and W. Metzner, Phys. Rev. 
Lett. 72, 316 (1994). 

LE. Dzyaloshinskii and A.l. Larkin, Zh. Eksp. Teor. Fiz. 

61, 791 (1971) [Sov. Phys. JETP 34, 422 (1972)]. 

A. I Larkin and Yu. N. Ovchinnikov, Zh. Eksp. Teor. 

Fiz. 55, 2262 (1968) [Sov. Phys. JETP 28, 1200 (1969)]; 

Nonequilibrium superconductivity, (Elsevier, Amsterdam, 

1984). 

'^^ G. Eilcnbcrgcr, Z. Phys. 214, 195 (1968). 

^° J. Schwinger, Phys. Rev. 128, 2425 (1962). 

A.O. Gogolin, A.A. Nersesyan, and A.M. Tsvelik, 
Bosonization and Strongly Correlated Systems, (Cam- 
bridge University Press, UK, 1998), Chap. 14. 
K.B. Efetov, Phys. Rev. Lett. 79, 491 (1997); Phys. Rev. 
B 56, 9630 (1997). 

L.P. Gorkov, Zh. Eksp. Teor. Fiz. 34, 735 (1958) [Sov. 
Phys. JETP 7, 505 (1958)]; Y. Nambu, Phys. Rev. 117, 
648 (1960). 

®* F. Essler and A. A. Nersesyan, Private communications. 

5. Lukyanov, Nucl. Phys. B 522, 533 (1998). 

^® N. Hatano and D. Nelson, Phys. Rev. Lett. 77, 570 (1996) 
^'^ K.B. Efetov and l.L. Alcincr (work in progress). 
'^^ K.G. Wilson and J. Kogut, Phys. Rep. C12, 75 (1974). 
A.A. Abrikosov and I.M. Khalatnikov, Sov. Phys. JETP, 

6, 888 (1958). 

®° W. Kohn and J.M. Luttingcr, Phys. Rev. Lett., 15, 524 
(1965); J.M. Luttinger, Phys. Rev., 150. 202 (1966). 



